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ABSTRACT 

A  determination  was  made  of  the  vapor  pressure  and  heat  of  vaporization 
of  solid  nickel  between  1233  and  1658°K.  The  value  of  the  second  law 
heat  of  vaporization  which  was  obtained  was:  (AHg)y  =  115.8  ±5.5 
kcal/mole,  while  the  value  of  the  third  law  heat  was:  (&Hg)v  =  102.6 
±2.1  kcal/mole.  A  statistical  treatment  of  all  the  available  data  of 
the  vapor  pressure  of  nickel  gave  the  following  "best  value"  for  tla 
heat  of  vaporization  of  solid  nickel-  (aHq)v$  =  101.45  ±0.62  kcal/irole. 

The  best  vapor  pressure  equation  of  solid  nickel  between  1250°K  and 
1710°K  which  was  obtained  was:  In  Ps(mm)  =  23.50  ±0.35  -  -—--y 
while  the  vapor  pressure  equation  for  liquid  nickel  between  1816  and 
1895°K  was:  In  P^(mm)  =  21.98  ±0.05  -  ±2-^.  Corrected  estimates 

of  the  heats  of  dissociation  of  the  molecules  (Nig)  and  (NiO)  were  also 
obtained  and  were:  Dg(Nig)  =  61.7  ±0.6  kcal/mole;  and  Dg(NiO)  =  92.7 
±0.9  kcal/mole  (millimeter  equals  Torr). 


Best  values  of  the  absolute  ionization  cross  sections  to  be  used  in 

mass  spectrometry  for  Ag,  0,  0g,  and  Ni  were  deduced  from  the  literature. 

The  maximum  cross  section  for  silver  at  72  volts:  cm72  (Ag)  =  4  63  x 

10-16  cm2;  while  o5Q  (Ag)  =  4.48  x  10"16  cm2;  ag0  (Ag)  =  4.50  x  10"lb 

cm2;  the  maximum  cross  section  for  0  atom  at  87.8  volts:  cr^gg^  = 

1.56  x  10“16  cm2;  the  maximum  (total  cross  section)  for  0?  at 

-1  fi  ?  ^ 

123  volts:  (amlgg)  =  ^-80  x  ^  cm'  .  The  maximum  cross  section 

t  -16  2 

for  nickel  at  approximately  70  volts  is:  ^^(Ni'  -  4.19  x  10”  cm 

v’hile  the  cross  section  at  60  volts  is:  agg(Ni)  =  4.10  x  10“16  cm2. 
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INTRODUCTION 


Since  there  is  a  continuous  need  for  information  on  the  thermodynamics, 
chemical  kinetics,  and  the  physical  behavior  of  materials  at  high 
temperatures ,  studies  of  the  mass  loss  rates  and  kinetics  of  phase 
transformations,  including  the  kinetics  of  vaporization  of  some 
selected  materials  have  been  undertaken.  To  increase  our  capability 
to  carry  out  these  studies  the  temperature  range  of  the  experimental 
techniques  which  were  used  in  previous  mass  loss  studies  has  been 
extended.  The  system,  which  was  constructed  to  dete*  nine  mass  loss 
rates,  was  characterized  by  establishing  the  operati ng- power  levels  of 
the  induction  furnace,  and  calibrating  the  techniques  for  measuring 
temperatures;  and  also  by  determining  the  reproducibility  and  errors 
of  mass  loss  measurements.  The  material  chosen  for  the  characterization 
studies  was  nickel. 

Upon  examination  of  available  information  on  the  heats  of 
vaporization  and  the  vapor  pressure  of  nickel,  it  became  apparent  that 
a  need  existed  for  a  determination  of  the  best  value  for  these  quantities; 
therefore,  a  statistical  evaluation  of  the  information  in  the  literature 
was  undertaken.  A  recalculation  of  the  values  of  the  vapor  pressures 
obtained  by  the  mass  spectrometric  technique  was  also  performed  due  to 
availability  of  more  recent  data  on  ionization  cross  sections  and 
thermodynami c  fimcti ons . 


The  report  is  divided  into  sections  describing  the  experimental 
technique  used  in  this  work,  and  the  resulting  measured  vapor  pressure 
and  temperature  relation  of  nickel;  the  statistical  evaluation  of  the 
data  available  in  the  literature  on  the  vapor  pressures  of  nickel  and 
heats  of  vaporization  of  nickel;  the  reestimation  of  the  heats  of 
vaporization  of  nickel  oxide;  reestimation  of  the  heats  of  dissociation 
of  the  NiO  and  Ni 2  molecules  and  a  determination  of  the  best  available 
values  for  the  ionization  cross  sections  of  silver  and  nickel  atoms. 


Throughout  this  report  mm  and  Torr  are  used  interchangeably. 
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PART  I 

EXPERIMENTAL  PROCEDURES 
A.  Description  of  the  Apparatus 

To  determine  vapor  pressures  of  materials  at  high  temperatures  it 
is  desirable  to  make  continuous  weighings  of  the  samples  without  exposing 
them  intermittently  to  the  atmosphere  since  exposure  to  the  atmosphere 
often  results  in  absorption  of  gases  which  cause  errors  in  the  determination 
of  the  mass  loss  rates.  The  errors  from  this  source  were  reduced  by  the 
use  of  a  recording  electrobalance,  coupled  with  an  induction  furnace, 
and  heating  the  sample  in  vacuo.  A  sample  thus  heated  can  be  continuously 
or  intermittently  weighed.  To  utilize  the  apparatus  properly, 
observations  were  made  on  the  temperature  uniformity  of  the  heated 
sample,  on  the  variations  of  the  recorded  weight  of  the  sample  with  the 
power  level  of  the  furnace,  and  on  the  sensitivity  of  the  microbalance. 

The  system  was  tested  by  using  it  to  measure  the  vapor  pressure  of  nickel 
by  determining  the  mass  loss  rates  of  a  cylinder  of  nickel. 

The  apparatus  consisted  of  a  6  in  i.d.  x  11  in  pyrex  bell  jar 
(a)  Figure  1  (photographs.  Figures  2  and  3)  which  housed  the  RF 
induction  coil  (b)  and  specimen  (c),  which  was  mounted  directly  over 
the  diffusion  pump  (h)  to  afford  maximum  pumping  speed  from  the 
specimen  chamber.  The  coil  (b)  consisted  of  five  turns  of  1/4 -inch 

copper  refrigeration  tubing  and  was  2-1/4  inches  i.d.  x  1-3/4  inches  in 
length.  The  specimen  (c)  was  suspended  and  centered  in  the  inside  of 
the  coil  (b)  by  means  of  a  0.005  inch  diameter  tungsten  wire  looped 
around  the  specimen  so  as  to  form  a  harness;  the  two  ends  were  hooked  to 
an  inconel  wire  three  inches  above  the  specimen,  and  the  inconel  wire, 
in  turn,  was  hooked  to  a  nylon  thread  hanging  from  the  end  of  the  balance 
beam  at  { d) .  The  electrobalance  (e)  (a  Cahn  type  RH)  was  enclosed  in 
a  glass  vacuum  bottle  (f).  The  balance  was  coupled  to  the  bell  jar  and 
pumping  system  via  a  short  length  of  vacuum  hose  (g)  which  afforded  a 
flexible  connection.  The  hose  was  raised  to  allow  insertion  or  removal 
of  the  specimen  from  the  system.  The  entire  system  was  evacuated  to 
Itf5  mm  pressure  by  means  of  a  6-inch  oil  diffusion  pump  (h)  and  a  Welch 
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Figure  1.  Experimental  Apparatus  for  Measuring  Weight  Losses 
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Figure  2.  Heating  Coil  and  Micro  Balance  as  Noted  in  Figure  1 
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Model  1397  Duo  Seal  mechanical  vacuum  pump  (i).  Baffles  were  inserted 
as  shown  at  (j)  in  order  to  reduce  the  back  diffusion  of  the  pump  oil 
into  the  specimen  chamber.  The  brass  flange  (k)  was  hollow  and  water 
was  pumped  through  it  for  cooling.  The  outside  of  the  bell  jar  (a) 
was  surrounded  by  a  water  jacket  through  which  was  circulated  distilled 
water  to  cool  the  specimen  chamber  surfaces  and  thus  remove  condensable 
gases.  The  cooling  water  was  circulated  in  a  closed  system  to  avoid 
corrosion  of  the  equipment,  and  the  formation  of  deposits  on  the  bell 
jar  surfaces.  The  induction  furnace  (1)  and  RF  coil  (b)  were  cooled  by 
the  same  system.  A  Leeds  and  Northrup  disappearing  filament  type 
optical  pyrometer  (m)  was  used  to  read  specimen  temperatures  via  the 
mirror  (n). 

The  induction  furnace  (1)  was  a  Sealomatic  Electronics  Corp.  Model 
2000-20KW  which  was  operated  at  '3  megaherz. 

The  electrobalance  (e)  was  monitored  by  means  of  a  Hewlett  Packard 
Model  419A  D.C.  null  voltmeter  (o).  The  voltmeter  received  its  input 
from  the  balance  control  box  (p).  The  output  of  the  voltmeter  (o)  was 
fed  to  a  Brown  (Model  153)  single  point  recorder  (q)  when  continuous 
weight  change  determinations  were  performed.  The  pressure  of  the  system 
was  determined  by  means  of  an  ionization  gauge  (r).  The  gauge  (r)  and 
pressure  control -monitor  box  (s)  were  Veeco  Vacuum  Gauqe  type  RG  31x. 

The  power  input  to  the  specimen  was  manually  adjusted  at  furnace  control 
knob  (t).  Power  input  and  stability  were  monitored  at  the  oscillator 
tube  plate  current  and  voltage  meters.  The  electrobalance  was  mounted 
on  a  3/4-inch  thick  plywood  platform  contained  in  an  aluminum  framework. 
One  end  of  the  framework  was  anchored  to  a  concrete  sidewall  of  the 
building  while  the  other  end  was  suspended  from  a  concrete  ceiling  by 
means  of  two  1/2-inch  steel  tie  rods  (see  photograph  Figure  3).  The 
mounting  minimized  vibrations  and  isolated  the  balance  from  the  induction 
furnace. 
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B.  Experimental  Procedure 

(1)  Calibration  of  Pyrometer  with  NBSa-STD  lamp  (IPTS  1948) 

The  optical  pyrometer  was  calibrated  against  a  standard  ribbon 
filament  strip  lamp  at  seyeral  temperatures  between  900  and  1600°C. 

The  temperatures  read  or-  the  pyrometer  during  cal i brat were  found  to 
be  within  ±5°C  of  the  standard  lamp  certification  curve.  The  pyrometer 
was  also  "calibrated"  by  placing  the  standard  lamp  in  "situ"  of  the 
sample;  thus  correcting  for  radiation  losses  due  ta  absorption  and 
length  of  path  between  sample  and  pyrometer.  These  calibrations  were 
carried  out  by  placing  the  standard  lamp  at  position  (b)  (Figures  1  and 
4)  and  focusing  the  pyrometer  on  the  image  of  the  lamp  in  mirror  (n). 
The  pyrometer  was  positioned  at  (m).  Two  observers  made  sets  of  three 
temperature  readings  each  on  the  pyrometer  for  each  fixed  power  level  of 
the  lamp,  resulting  in  six  temperature  readings  for  each  calibration 
point.  The  two  sets  of  readings  at  any  given  temperature  were  foun,.  to 
agree  to  ±5°C.  me  calibration  curve  of  "true"  vs.  "read"  temperatui  ? 
consisted  of  six  different  points  thus  determined. 


14" 


je"  dlo.  x  deep 
black  body  cavity 


Cylindrical  Specimen 


Figure  4.  Arrangement  of  Optical  Pyrometer  for  Temperature 
Measurements 
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(2)  Weighings 

The  masses  of  the  samples  were  determined  by  reading  the  mass 
dial  of  the  microbalance  control  box  (d)  when  voltmeter  (o)  indicated 
zero  voltage,  each  division  on  the  mass  dial  corresponded  to  0.5  milligram 
and  the  dial  could  be  read  to  two-tenths  of  one  division  giving  a 
sensitivity  of  0.1  milligram.  The  reproducibility  of  individual 
weighings,  however,  was  found  to  be  ±0.3  milligram,  on  the  basis  of 
thirteen  individual  weighings  of  a  specimen  under  experimental  conditions. 
The  electrical  noise  level  was  low,  and  no  appreciable  fluctuations  of 
the  meter  (o)  were  observed  when  the  mass  dial  was  read.  The  electro¬ 
balance  was  checked  prior  to  tech  experiment  by  the  addition  and 
subtraction  of  milligram  calibre tion  weights  to  the  balance  tare  weight. 

The  weight  as  read  out  0'<  the  mass  dial  it  (p)  agreed  with  the  calibration 
weight  to  ±0.1  milligram. 

(3)  Specimen 

The  nickel  specimen  was  machined  from  a  one-half-inch  diameter 
nickel  rod  of  99.95%+  purity  which  was  obtained  from  Galland  Schlesinger 
Chemical  Company,  Carle  Place,  New  York.  A  black  body  cavity  1/16  inch 
in  diameter  by  1 /4-inch  deep  was  drilled  'n  the  center  of  one  end  of  the 
specimen  such  that  it  could  be  sighted  for  the  optical  temperature 
measurement  (see  Figure  4).  The  specimen  was  cleaned  with  benzene  to 
remove  any  oil  or  other  impurities  which  might  be  absorbed  on  the  surface, 
and  the  cylinder  was  suspended  as  described  in  Section  A.  The  required 
counterweights  were  added  to  pan  (v)  so  that  the  total  weight  loss  could 
be  read  out  at  (p)  without  breaking  the  vacuum.  Next,  the  specimen 
was  centered  as  well  as  possible  in  coil  (b)  and  the  system  evacuated  to 

C 

10”  Torr  and  power  applied  from  induction  heater  (1)  by  adjusting  the 
control  (t)  to  the  desired  power  level.  The  sample  was  outgassed  by 
heating  to  approximately  1400°C  for  15  minutes  prior  to  each  run.  The 
power  output  as  observed  at  the  meters  (u)  was  constant  to  ±4%  at  all 
temperatures  throughout  the  runs.  After  the  sample  reached  the  desired 
temperature  the  sample  weight  was  observed  at  the  recorder  (q),  until  a 
constant  weight  loss  rate  was  observed  indicating  that  the  sample  was  at 
a  constant  temperature  and  had  been  adequately  degassed.  The  initial 
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time  of  a  run  was  taken  as  the  time  at  which  the  mass  loss  rate  became 
constant. 

When  the  power  was  turned  on,  the  nickel  specimen  experienced  a 
solenoid  force  which  was  due  to  the  field  in  the  induction  coil.  This 
force  varied  with  the  power  setting,  but  was  constant  at  a  given  powjr 
setting  as  is  indicated  by  the  agreement  in  waight  loss  rates  as 
determined  from  the  ,otal  weight  loss  measurements  divided  by  the  time 
and  ‘•'lat  determined  from  the  recorder  chart  slope  of  Wt  *ght  vs  time 
measurements.  The  agr*a:ent  was  within  12?  when  nit/  terminations 
were  compared.  The  solenoid  force  was  of  the  order  ».  _ne  milligram  ?t 
the  average  power  setting  of  about  0.3  kilowatt. 

The  sample  was  heated  at  a  given  temperature  long  enough  so  that  the 

resulting  weight  loss  was  several  milligrams.  Thus  keeping  the  probable 

weighing  error  to  less  than  1?  without  causing  an  appreciable  r.nanue  in 

2 

the  sample  surface  area,  the  sample  area  was  between  5  and  6  cm  .  The 
vaight  of  the  specimen  was  recorded  immediately  before  and  after  each 
heating  period  at  a  given  temperature  with  the  heater  powe.  turned  off; 
thus  a  total  v.eight  loss  at  each  temperature  was  obtained. 

Weight  losses  were  also  determined  from  a  continuous  weighing  of  the 
sample  as  a  function  of  time  when  the  power  settings  were  0.7  kw  or  le<s 
at  a  frequency  of  2.8  Megaherz.  At  power  settings  higher  than  0.7  kw 
the  noise  level  was  too  high  to  give  a  consistent  weight  with  th»* 
on  and,  therefore,  the  only  weight  determinations  which  were  made  were 
those  obtained  with  the  power  shut  off. 

For  a  given  determination  of  the  heat  of  vaporization,  mass  loss 
rates  wera  measured  at  six  to  ten  temperatures  spaced  at  equal  intervals. 
To  obtain  adequate  mass  losses  at  low  temperatures,  runs  lasting  several 
hours  were  necessary,  whereas,  at  the  high  temperatures,  runs  of  twenty 
minutes  gave  adequate  mass  losses.  All  runs  were  at  least  twenty 
minutes  long  in  order  to  make  the  error  in  the  time  determination 
insignificant,  and  to  minimize  the  error  in  mass  loss  determinations 
which  occur  during  the  time  period  that  the  sample  is  being  heated  and 
cooled. 
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C.  Treatment  of  Data 


If  the  rate  of  mass  loss  per  unit  area  from  a  freely  evaporating 
solid  is  equal  to  the  .iapori7at’.on  rate  of  the  solid  in  eqai librium  with 
its  vapor,  then  rate  of  mass  loss  per  unit  area  (m-)  due  to  a  giver, 
gaseous  species  (i)  at  the  temperature  (T)  ic  related  to  the  standard 
heat  of  vaporization  [AH°-]  for  the  soecies  (i)  by  Equation  1  (Reference  1). 

d(lnmj)  d(InFj)  T  AH°vi  T 

*d(i7r)  "  d(l/T)  *  2  R  *  2  (1) 

In  Equation  1,  P.  is  the  equilibrium  vapor  pressure  of  the  ith  species. 

The  slope  of  a  plot  of  In  m.  vs.  1/T,  is  proportional  to  the  neat  of 
vaporization  if  we  /'o-'ime  tnat  the  difference  in  the  neat  capacity  of 
the  gas  and  solid  is  zero.  The  error  due  to  the  differences  in  the  heat 
capacity  is  minimized  if  we  assume  that  AH°  is  ti  ken  at  the  mean 
temperature  of  the  determi nation. 

In  considering  the  vaporization  of  nickel,  the  gaseous  species  is 
assumed  to  be  only  nickel  atoms;  (References  2,  4).  Therefore,  the 
vapor  pressure  can  be  calculated  by  using  the  following  relation: 


•  -2  -1 
if  m  =  g  x  cm  x  sec 

P  =  pressure  in  millimeters 
M  =  molecular  weight  gm/mol 
T  =  degrees  kelvin  for  nickel,  M  =  58.71  g/mol 

and  P  = 


m  yr 

0.447 
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A  computer  program  'as  written  to  obtain  a  least  square  fit  to  a 
linear  function  of  Tc.  m  vs.  1/T.  The  data  used  in  the  program  was  the 
observed  time,  the  mass  loss,  and  the  temperature.  The  slope 
[din  m/d  (1/T)3»  F  ;ure  5.  is  related  to  [AH® J  by  Equation  1  above. 

The  value  of  [AH°j  f/>rs  obtained  was  referred  to  ?,  the  temperature 
obtained  by  averaging  th  .*  Averse  temperatures  of  the  set  of  observations 
of  the  mass  loss  vs.  time,  and  is  designated  in  the  literature  as  the 
2nd  law  heat  of  vaporization.  The  series  of  values  for  the  second  law 
heat  of  vaporization,  f.AH;VJr»  thus  obtained  were  corrected  to  the 
reference  stater,  at  the  reference  temperature  [T  ]  of  298°K  and  0°K  by 
the  use  of  Equation  3  below.  (Reference  3). 

[H-T-H;r]s-[HnL"NTsHr  is) 

The  best  values  of  AHg  and  AH|gg  which  were  obtained  are  shown  in  Table 
II  and  Figure  5. 

Third  law  values  of  £AH°]j  were  obtained  utilizing  free  energy 

r 

functions  as  given  in  the  standard  thermodynamic  tables  (Reference  3). 

(G°  -  HSgJ 

For  example  values  of - =p£2-  f  r  the  gas  and  solid  were  used  as 

shown  below,  in  Equation  5  to  obtain  the  [AH°]  for  the  vaporization  of 
the  solid  nickel.  Reaction  4, 

(Ni)$=  (Ni)y  (4) 


j 

l 

1 

Table  II  gives  the  data  used  to  obtain  values  of  the  third  law  heat  of  * 

valorization. 
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An  estimate  v*as  made  of  the  probable  errors  in  each  of  the 
determinations  of  the  vapor  pressures  on  a  basis  of  the  provable  errors 
in  the  measurements  (mass  loss,  and  temperature.  Table  I).  The 
estimated  ,'.r*'or  in  the  weighing  was  ±0.3  milligram  thus  the  %  error 
varies  with  the  totzl  mass  loss  (Am  Table  IV)  at  a  given  point.  This 
error  was  compared  v/ith  that  determined  from  a  least  squares  computer 
program  which  determined  the  error  in  the  slope  of  the  line*  In  P  vs. 
1/T,  as  shown  in  Figure  r,.  The  error  in  the  slope  can  be  compared  to 
those  calculated  on  the  basis  of  the  probable  errors  in  T  and  Am,  by 
comparing  the  errors  in  the  pressure  obtained  by  the  equations  given  in 
Table  IV.  Comparison  of  errors  was  made  using  the  data  in  Table  II. 

The  results  are  given  in  Table  IV. 

The  data  reported  on  the  vapor  pressures  of  nickel  and  heats  of 
vaporization  of  nickel  as  determined  by  the  second  and  third  law  methods 
were  combined  with  other  data  in  the  literature  to  determine  the  best 
values  for  the  vapor  pressures  and  heats  of  vaporizati'”’.  Kart  II, 

Sec.  1  (3(1),  of  this  report  discusses  the  statistical  methods  for 
combining  data  from  various  laboratories  to  obtain  the  best  values. 
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TABLE  I 


MASS  LOSS  RATES  FOR  LANGMUIR  VAPORIZATION  OF  NICKEL 


At 

Am 

m  6 

Pressure 

Run  No. 

T°K 

min. 

m<j 

g/cnr^/sec  x  10 

Torr  x  10° 

3a,b.c 

1233 

126 

0.1 

0.0033 

0.2598 

4b, c 

1658 

8 

50.6 

27.00 

2443. 

4 

1563 

10 

11.6 

4.03 

356.7 

4 

1393 

102 

2.0 

0.73 

6.083 

6d 

1558 

29 

22.3 

2.143 

189.2 

5 

1447 

59 

9.5 

0.4487 

38.18 

6 

1513 

26 

31.6 

3.388 

294.8 

6 

1277 

248 

0.45 

0.0051 

0.4078 

6 

1345 

169 

0.55 

0.0091 

0.7465 

6 

1428 

117 

1.05 

0.0253 

2.139 

6 

1546 

43 

10.9 

0.7148 

62.83 

6 

1648 

11 

26.7 

6.846 

621.8 

6 

1562 

34 

16.8 

1.393 

123.2 

6 

1610 

21 

3~.0 

£.700 

421.8 

lOe 

7  <-58 

142 

5.3 

0.118 

10.08 

10 

1488 

85 

6.3 

0.235 

20.28 

10 

1523 

44 

6.8 

0.490 

42.80 

10 

1563 

30 

9.2 

0.970 

85.79 

10 

1593 

18 

17.7 

3.11 

277.6 

10 

1613 

12 

14.5 

3.82 

343.1 

10 

1613 

20 

22-3 

3.53 

317.0 

lie 

1403 

339 

2.5 

0.0237 

1.986 

11 

1373 

394 

1.7 

0.0139 

0.5150 

11 

1363 

274 

0.9 

0.0106 

0.8760 

12e 

1353 

215 

0.3 

0.0045 

0.3310 

12a, b 

1483 

48 

4.5 

0.302 

26.00 

12a, b 

1543 

15 

4.8 

1.03 

90.51 

12 

1433 

59 

0.5 

0.027 

2.270 

a  These  points  eliminated  from  2nd  law  determ1  nation, 
b  These  points  eliminated  from  3rd  law  determination, 
c  Specimen  surface  area  =  4.05  cm2, 
d  New  specimen  same  stock  as  c.  Area  =  5.91  cm. 
e  New  specimen  same  stock  as  c.  Area  =  5.17  cm2. 
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TABLE  II 

THIRD  LAW  VALUES  (AHg)y  FOR  THE  VAPORIZATION  OF  NICKEL 


T°K 

Pressure 
atmospheres 
x  109 

AH°298 

kcal/mol 

1345 

0.982 

102.92 

1428 

2.810 

106.16 

1546 

82.70 

97 . 75a 

1447 

50.20 

99.25 

1558 

249.0 

101.72 

1562 

162.0 

103.33 

1563 

469.0 

100.06 

1610 

555  0 

102.46 

1648 

818.0 

104.17 

1658 

3210.0 

99.65 

1363 

1.153 

103.86 

1353 

0.435 

105.74 

1433 

2.990 

106.37 

1483 

34.20 

102.83 

1543 

119.1 

102.46 

1523 

56.30 

104.00 

1488 

26.70 

103.86 

1458 

13.30 

103.84 

1613 

417.0 

103.57 

1613 

dR1.0 

103.34 

1593 

066.O 

102.74 

1563 

112.=# 

104.50 

1403 

2.620 

105.17 

1373 

0.6780 

106.03 

1393 

2.004 

100.72b 

1233 

0.3413 

97.09ab 

'.211 

0.5366 

99.36b 

1513 

387.9 

97 . 53ab 

a.  Points  not  used  to  determine  3rd  law  value. 

b.  Points  not  used  to  determine  2nd  law  value. 

3rd  law  AHq  =  102.6  ±2.1  kcal/mole. 

3rd  law  AHggg  =  103.1  ±2.1  kcal/mole. 

2nd  law  AHq  =  115.8  ±5.5  kcal/mole. 

2nd  law  AHggg  =  116.3  ±5.5  kcal/mole. 
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TABLE  III 

FREE  ENERGY  FUNCTIONS  FOR  NICKEL  SOLID  AND  NICKEL  GAS  (a) 


[-(S°-H"st)/T]s 

C-(G°-H°st)/T]( 

T°K 

cal/°K/ir.ole 

cal/°K/mole 

1345 

47.74 

12.43 

1428 

48.02 

12.80 

1546 

48.38 

13.66 

1447 

48.08 

12.88 

1558 

48.43 

13.35 

1562 

48.44 

13.36 

1553 

48.44 

13.37 

1610 

48.58 

13.56 

1648 

48.69 

13.71 

1658 

48.72 

13.75 

1363 

47.81 

12.51 

1353 

47.78 

12.46 

1433 

48.04 

12,81 

1483 

48.20 

13.02 

1543 

48.38 

13.66 

1523 

48.32 

13.20 

1488 

48.21 

13.06 

1458 

48.12 

12.93 

1613 

48.59 

13.57 

1613 

48.59 

13.57 

1593 

48.53 

13.49 

1563 

48.44 

13.37 

1403 

47.94 

12.24 

1373 

47.84 

12.56 

1393 

47.91 

12.65 

1233 

47.48 

11.92 

1277 

47.51 

12.12 

1513 

48.29 

13.16 

*st  -  standard  state  at  298°K. 

a.  Tabulated  values  interpolated  from  Hultgreen}  Reference  3 
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TABLE  IV 

ESTIKATED  EXPERIKEKTAL  ERRORS  IN  THE  OBSERVED  VAPOR  PRESSURES  OF  NICKEL, 
AND  DEVIATIONS  BETWEEN  OBSERVED  AND  CALCULATED  PRESSURES 


%  Error 

Observed 

Predicted 

Calculated  . 

P(Torr)  x  10b  D 

%  dev. 

T°K 

P(Torr}  x  10° 

in  P(a) 

1345 

0.7465 

33.6 

0.1493 

+400.0 

1353 

0.3310 

100.0 

0.7580 

+56.3 

1363 

0.8760 

33.0 

0.4380 

+100.0 

1373 

0.5150 

17.7 

1.772 

-70.9 

1403 

1.990 

12.0 

2.388 

-16.7 

1428 

2.139 

29.5 

11.29 

-81.1 

1433 

2.270 

60.0 

14.44 

-84.3 

1447 

38.18 

3.4 

5.83 

+555.0 

1458 

10.10 

5.7 

10.04 

+0.6 

1483 

26.00 

6.7 

13.47 

+93.0 

1488 

20.30 

4.8 

25.17 

-19.4 

1523 

42.80 

4.4 

73.19 

-30.4 

1543 

90.50 

6.3 

85.38 

+6.0 

1546 

62.88 

3.0 

161.4 

-61.0 

1558 

189.2 

1.5 

74.2 

+155.0 

1562 

123.2 

0.7 

158.9 

-22.5 

1563 

85.80 

3.3 

272.8 

-68.6 

1563 

356.7 

1.9 

42.7 

+735.4 

1593 

278.0 

1.7 

289.1 

-3.S 

1610 

421.8 

1.0 

404.4 

+4.3 

1613 

343.0 

2.1 

617.4 

-44.4 

1613 

317.0 

1.4 

684.7 

-53.7 

1648 

621.8 

1.3 

1623. 

-61.7 

1658 

2442. 

0.8 

442 

+453.0 

(a) 

Predicted  errors  on 
measurement  of  P. 

the  basis 

of  expected  errors 

in  the 

%  error  =  jjjjp  +  - 

x  100 

(b)  Calculated  on  basis  of  least  squares  fit  as  shown  in  Figure  5. 

(c)  %  deviation  from  least  squares  line. 

[observed  P  -  calculated  P  lnR 
calculated  P  x  iuu 
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in  P._  .  =  27.8$  - 

(Torr) 


5.740  x  10* 


end  line  Intercepts : 


in  P  =  -  0.837 


-20.00  I 


Figure  5.  Secom'  Law  Plot  for  the  Determination  of  (AHg)y  of 

Nickel  (solid),  (a)  Po’nts  with  deviations  >2a  (r.m.s) 
deviation  eliminated. 
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PART  II  - 

EVALUATION  OF  THERMODYNAMIC  DATA  ON  THE  VAPORIZATION 
NICKEL  AND  NICKEL  OXIDE 

SECTION  1  Examination  of  the  uata  on  the  Vapor  Pressures  and 

Heats  of  Vaporization  of  Nickel 

An  opportunity  to  determine  the  best  value  of  t'  e  he it  of  vaporization 
of  nickel  and  to  obtain  an  expression  for  the  vapor  pressure  l/as  offered 
by  the  available  data  in  the  literature.  The  variety  of  techniques 
reported  for  the  determination  of  the  vapor  pressure  of  nickel  all 
reduce  to  a  determination  of  the  mass  loss  rate  from  a  sample  as  a 
function  of  temperature,  with  the  exception  of  the  determination  by 
Grimley  et  al  (Reference  4).  All  the  data  can  be  reduced  in  terms  of 
■‘■be  second  and  third  laws  of  thermodynamics  as  discussed  under  "Treatment 
Data",  Part  I,  Section  1C  of  this  report. 

In  determining  the  suitability  of  data  for  use  in  determining  the 
best  value,  certain  criteria  were  applied  for  evaluating  the  experimental 
techniques  and  also  an  examination  was  made  of  the  internal  and  external 
consistencies  of  the  data. 

The  following  discussions  cover  each  reported  set  of  experiments 
and  also  the  opinions  of  the  validity  of  some  of  the  data  as  expressed 
by  other  evaluators.  The  sources  of  data  are  given  in  tiic  references 
and  in  Table  V. 

A.  Evaluation  of  Experimental  Techniques 

In  the  evaluation  of  the  experimental  techniques  we  took  note  of 
(1)  methods  of  temperature  measurement,  (2)  condition  of  samples, 

(3)  method  of  weight  loss  determination,  (4)  consistencies  between 
second  and  third  law  values  of  heats  of  vaporization,  (5)  external 
consistencies  of  heats  of  vaporization. 

a.  Haury,  this  report,  Part  I,  determined  the  temperatures  of 
his  samples  by  sighting  on  a  black  body  hole  in  the.  sample  with  a 
calibrated  pyrometer.  The  complete  experimental  calibration  of  the 
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TABLE  V 


REPORTED  VAPOR  PRESSURES  OF  NICKEL  USED  IN  DETERMINING 
BEST  VALUES  OF  BY  SECOND  AND  THIRD  LAM  METHODS 


Temperature 

°K 

Pressure 

Torr 

Temperature 

°K 

Pressure 

Torr 

Temperature  Pressure 

°K  Torr 

Haury, Table  I 

1393c,d 

Johnston  6  Marshall, 

Langmuir  £  McKay,  Reference 

0.6083*10  ° 

Reference  (8) 

1318 

0. 2436x1 0-1* 

13456 

0.7465x10-5 

1583 

0.3520x10-6 

1409 

0.5072x10-4 

1610 

0.4218x10-3 

1308 

0.4260x10-6 

1532 

0.3301x10*3 

1546 

0.5288x10-4 

1307 

0.4560x10-6 

1544 

0.3850x10-3 

1447 

1563 

0.3818x10-4 

0.3567x10-3 

1397 

1397° 

0.5500x10-5 

0.6200x10-5 

16046 

0.1271x10*2 

1558 

0.1892x10-3 

1387 

0.4330x10-5 

Nesmeynov  S  Teh  Tik-Hang, 

1562 

0.1232x10-3 

1415 

0.8280x10-5 

References  (5),  (6)e 

1658 

0.2443x10-2 

1465 

0.2830x10-4 

1525 

0.1170 

0.4620x10-' 

0.4690x10-1 

0.2330x10-1 

0.2040x10-1 

1648 

0.6218x10-3 

1507 

0.7540x10-4 

1483 

1277C 

0.4078x10-8 

157p 

0.3440x10-3 

1478 

1428 

0.2139x10-5 

1078 

0.3460x10-3 

1443 

1403 

U.  1990x10-5 

1397c 

0.6200x10-5 

1439 

1373 

1363 

0.5150x10-6 

0.8760x10-6 

1466c 

0.3290x10-4 

1401 

1411 

0.9800x10*; 

0.1060x10-' 

1433 

0.2270x10-5 

Grimley,  Burns  t.  Ingrahm, 

1365 

0.3680x10-2 

1523 

0.4280x10-4 

Reference  (4)° 

!320d 

0.1640x10-2 

1488 

0.2030x10-4 

1575C 

0.2060x10  i 

1550“ 

0.1730 

1458 

1513 

0.1010x10, 

0.3170x10”, 

1596c 

1606 

0.3400x10”, 

0.4130x10”, 

1371c 

0.3540 

1S13 

0.3430x10”, 

1624 

0.6460x10”, 

0.6460x10”, 

Morris,  Zellars,  Payne  £ 

1593 

0.2780x10”]} 

1625 

Kipp, 

Reference  (11)'  . 

0.3870x10”} 
0.3860x10'} 
0.3320x10"} 
0.1280x10"} 
0.3060x10'} 
0.2640x10  } 
0.2210x10”} 
0.2020x10”} 
0.1910x10'} 
0.1780x10”} 
0.1750x10'} 
0.1510x10”} 
0.2040x10'} 
0.2430x10” 
0.1360x10”} 
0.2560x10"' 

1563 

0.8580x10”’ 

1630 

0.7170x10”, 

1895 

1543 

0.9050x10”’ 

1646 

0.1060xl0”‘ 

1894 

1353 

0.3310x10”° 

1651 

o.ii60xio”; 

1885 

14M  . 

1233*;’ j 

0.2600x10"’ 

1657 

0.1280x10”; 

1816 

0.2598x10”, 

1659 

0.1300x10”; 

1876 

1513c,d 

0. 2948x1 0”J 

1673 

1673 

0.1660x10”; 

0.1710x10”; 

1868 

1853 

Bryce,  Reference  (12)  , 

1679 

0.1910x10”; 

1848 

1273 

0.2450x10”, 

1684 

0.2130x10"; 

1842 

1252 

0.9800x10"' 

1679 

0.1930x10”; 

1839 

1423 

0.1 400x1 Op 

1684 

0.2140x10”; 

1836 

1391 

0.3690x10”° 

1707 

0.3130x10”; 

1828 

1341 

0.1100x10'° 

1709. 

0.3280x10”; 

1848 

1324 

1299 

0.1050x10”° 

0.7510x10”° 

1587d 

0.3360X10”-5 

1860 
1821  . 
1862° 

Uoca  points  were  omitted  whose  deviations  were  greater  than  2o  (r.m.s.  deviation)  from  second 
law  and  thrid  law  best  values. 

a. .  All  points  reported  included. 

b.  This  data  has  been  corrected  for  errors  in  cross  sections.  See  p.  22. 

c.  Eliminated  from  second  law  individual  determination. 

d.  Eliminated  from  third  law  individual  determination. 

e.  Not  used  in  final  composite  determination.  See  p.20. 

f.  Liquid  phase  vaporization  corrected  for  heat  of  fusion.  See  p.  21,  26. 

Not  used  to  determine  best  value  of  A  by  second  or  third  law. 

g.  Not  used  in  composite  third  law  determinations. 
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pyrometer  whiu..  Is  described  on  page  7  included  corrections  for  the 
absorption  of  radiation  by  the  sample  viewing  windows.  Although 
observations  of  the  "brightness"  temperature  on  the  external  surface  of 
the  cylinder  were  made,  no  data  was  recorded  on  the  extent  of  the 
variation  of  the  temperature  on  the  surface.  The  weight  loss  rate  was 
obtained  by  continuous  weighing  with  a  microbalance. 

The  lack  of  consistency  between  the  second  and  third  law  values  of 
the  heat  of  vaporization  which  is  shown  in  Table  II,  arises  from  the  fact 
that  the  second  law  values  are  subject  to  larger  errors  than  third  law 
values,  as  indicated  by  the  error  analysis  in  Part  I,  Section  1C. 

b.  The  data  of  Nesmevanov  and  Tik-Mang  (References  5,  6,  and  7) 
were  evaluated  but  not  utilized  because  the  vapor  pressures  which  were 
determined  by  measuring  the  weight  loss  rate  from  a  Knudsen  cell  were 

3 

10  times  larger  than  those  previously  and  subsequently  reported. 

c.  Johnston  and  Marshall  (Reference  6)  determined  the  vapor  pressure 
of  nickel  by  a  free  vaporization  or  Langmuir  technique.  They  determined 
the  temperature  by  sighting  on  a  black  body  hole  in  the  sample  which  was 
doughnut  shaped.  They  made  brightness  intensity  corrections  for  the 
absorption  of  radiation  by  the  windows.  The  shape  of  the  sample 
introduced  an  error  due  to  condensation  of  material  on  the  inside  of  the 
doughnut  hole;  a  correction  was  made  for  this  error.  The  weight  loss 
was  determined  by  the  use  of  a  sensitive  balance.  Although  the  authors 
claim  an  error  of  ±5%  in  the  pressures  measured,  the  data  does  require 
some  interpretation  since  there  appears  to  be  an  uncertainty  in  the  area 
of  the  condensation  surface.  The  data  presented  involved  an  error  of 
±5%  in  the  mass  loss,  and  ±5°K  in  the  temperature. 

d.  Langmuir,  Jones  and  MacKay  (References  9  and  10)  reported  the 
measurement  of  the  vapor  pressure  by  a  free  evaporation  from  nickel 
filament.  The  temperature  was  measured  with  a  pyrometer  and  an 
emissivity  of  0.36  was  assumed  for  the  nickel  at  temperature.  This 
assumption  leads  to  an  error  in  the  determination  of  the  temperature. 
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There  is  another  error  due  to  the  determination  of  the  length  of  wire 
at  temperature.  Although  estimates  on  the  limits  of  these  errors  could 
be  determined,  the  data  herein  falls  out  of  the  statistical  limits 
which  v.'ere  set  for  use  of  the  data  in  obtaining  the  best  values  of  the 
heats  of  vaporization  and  vapor  pressures  by  the  third  law  method. 

e.  Morris  et  al  (Reference  11)  measured  the  vapor  pressures  of 
liquid  nickel  by  a  gas  transport  method  utilizing  an  inert  gas,  such  as 
helium,  as  the  transport  gas.  The  metal  vapor  was  condensed  onto  the 
inner  surface  of  the  gas  outlet  tube.  The  gas  flow  rates  were  determined 
by  a  wet  gas  meter.  The  temperature  was  controlled  to  ±5°K  and  was 
determined  by  means  of  an  optical  pyrometer.  Corrections  for  absorption 
of  the  radiation  by  the  windows  were  made.  A  colorimetric  method  was 
used  to  determine  the  weight  of  nickel  deposited  by  the  gas  stream. 

The  possible  errors  in  the  spectroscopic  determination  of  the  nickel 
carried  by  the  gas  were  not  evaluated..  The  error  in  the  volume  of  the 
gas  gives  a  proportional  percentage  error  in  the  vapor  pressure.  The 
probable  error  in  the  gas  volume  is  less  than  1%. 

f.  Bryce  (Reference  12)  also  determined  the  vapo”  pressure  of  nickel 
by  the  vaporization  of  a  filament  using  a  target  to  catch  the  nickel 
vapor.  The  nickel  condensed  on  the  target  was  determined  o>  a  chemical 
technique.  No  further  data  is  available  on  this  technique. 

g.  Grimley  et  al  (Reference  41)  determined  the  vapor  pressure  of 
nickel  during  the. course  of  a  mass  spectrometric  investigation  of  the 
vaporization  behavior  of  nickel  oxide  by  combining  the  Knudsen  cell 
effusion  technique  with  a  mass  spectrometer.  The  temperature  was 
determined  by  sighting  on  black  body  holes  in  the  Knudsen  cell.  No 
information  was  given  on  the  possible  temperature  gradients  in  the  cell. 
There  was  some  evidence  that  a  reaction  took  place  in  the  Knudsen  cell 
which  was  made  from  aluminum  oxide.  There  is  no  evidence  that  this 
reaction  affected  the  measured  vapor  pressure  or  the  observed  heat  of 
vaporization  of  nickel. 
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Sines  a  mass  spectrometer  was  used  to  determine  the  vapor  pressure, 
a  discussion  of  this  technique  is  given  in  Appendix  B,  Section  2. 

Grimley  et  al  (Reference  4)  used  incorrect  ionization  cross  sections 
for  the  ionization  of  the  molecular  and  atomic  species  involved.  The 
error  in  the  ionization  cross  section  of  silver  as  noted  in  Appendix  B, 
Section  2,  introduces  an  error  in  the  geometric  factor  defining  the 
number  of  atoms  which  are  in  the  volume  defined  by  the  electron  beam. 

Other  errors  were  made  by  these  authors  in  reducing,  thei'  data  by  the 
third  law  technique,  for  example,  the  use  of  values  of  the  free  energy 
functions  of  NiO  which  were  inconsistent  with  available  molecular  data. 

A  discussion  of  this  paper  is  given  in  Part  II,  Section  2B. 

B.  Estimate  of  Second  and  Third  Law  Heats  of  Vaporization 

Although  a  statistical  technique  has  been  reported  (References  13 
and  14)  for  obtaining  "best  values"  of  "second"  and  "third  law"  heats  of 
vaporization  f»*om  vapor  pressure  data  obtained  from  a  series  of  experiments, 
and  from  a  numL »r  of  different  laboratories  or  sources,  in  this  report 
a  technique  (Reference  15)  proposed  by  Paule  and  Handel  was  utilized 
which  differs  noticeably  from  the  one  they  used  in  (References  13  and 
14). 


The  first  step  in  obtaining  the  "best  value"  was  to  examine  the  data 
of  the  individual  laboratories  for  internal  consistency,  and  also  for 
consistency  with  other  available  experimental  data,  and,  for  information 
about  the  experimental  techniques  used  as  noted  above,  p.  21.  In  the 
"second  law"  treatments  of  the  experimental  data,  the  data  from  each 
laboratory  was  fitted  by  a  least  squares  method  to  the  Equation  6 


In  P  -  A  +  — 
T 


(6) 


where  P  is  the  pressure  (Torr)  and  B  =  (AH|/R) ;  AH|  is  the  heat  of 
vaporization  at  the  temperature  T,  where  T  is  the  mean  inverse  temperature 

-  A-  T— "l  •  ’ 

calculated  by  T  =\ni  TJ  ;  T.  is  the  temperature  of  the  j-th  data 
point,  and  n.  is  the  number  of  data  points  reported  by  the  i-th  laboratory. 
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All  reported  pressures  having  a  deviation  from  the  calculated  values  of 
the  pressure  which  were  greater  than  2a  (a  =  r.m.s  'Jeviation)  v/ere 
eliminated  from  the  original  reported  data  in  determining  the  least 
squares  value  of  the  slope. 

(1)  Determination  of  Second  Law  (AHg)v,  8  and  A 

The  composite  best  value  of  (AHg)y  determined  by  means  of  the 
"second  law",  that  is,  the  slopes  o.f  the  various  curves  was  obtained 
by  first  correcting  the  determined  value  of  (AH|)  to  AHjj.  The  corrections 
were  made  by  using  the  enthalpy  tables  of  Stull  and  Sinke  (Reference  16), 
which  are  almost  identical  to  those  in  References  3,  17,  and  18,  since 
all  tables  round  o'f  to  the  same  values  to  the  second  decimal  place. 

The  list  of  data  used  in  this  determination  is  given  in  Table  V. 

Before  applying  the  "least  squares"  treatment  to  the  data  of  the 
individual  laboratories,  the  data  was  examined  for  signs  of  "systematic 
errors"  or  "non-random  errors"  by  arranging  the  deviations  of  the 
pressures  in  ascending  order  and  plotting  the  deviations  cn  probability 
paper.  In  all  cases  the  data  formed  several  groupings  which  could  be 
fitted  with  straight  lines.  The  deviations  of  the  plots  from  a  single 
straight  line  was  not  limited  to  the  end  points,  thus  indicating  that 
Ihe  errors  were  "non-random"  (Reference  19  and  20). 


The  best  values  of  A  and  B  of  equation  6  for  each  individual  set  of 
data  were  determined  by  successively  eliminating  data  points  whose 
deviations  between  the  measured  and  calculated  values  of  the  pressure 
were  greater  than  2 a  (r.m.s.  deviation).  In  some  cases  more  than  one 
recalculation  was  necessary,  thus  indicating  a  lack  of  statistical 
homogeneity  in  the  data. 


The  composite  best  values  of  A  and  B  were  obtained  by  utilizing  the 
statistical  treatment  described  in  (Reference  15)  to  determine  the 
applicable  weighing  factors  in  averaging  the  values  of  AHg  derived  by 
the  least  squares  treatment.  The  treatment  in  (Reference  15)  is  an 
attempt  to  weigh  the  effect  of  deviations  between  the  individual 
laboratory's  "best  values and  the  "within  the  laboratory"  deviations; 


the  individual  weighing 


factors  were  caiculav.ed 
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1.  To  calculate  the  quantity  defined  In  (Reference  15)  as  the  mean 
square  deviation  within  laboratory  data,  it  is  necessary  to  consider  the 
fact  that  r.m.s.  value  of  the  slope  is  determined  by  (n.-l)  quantities 
where  n-  is  the  number  of  data  points  reported  by  laboratory  i.  The 
mean  square  deviation  of  the  "within”  laboratories  d?-ta  is  defined  by 
the  quantity  MSy. 


MSW  = 


fin  ~2)erf 
N-k 


(7) 


2 

where  k  is  the  number  c*  laboratories,  o.  is  the  variance  of  the  slope 
calculated  for  laboratory  i,  (n.-2)  is  the  number  of  degrees  of  freedom 
in  determining  the  variance  of  the  slope,  and  N  =  ^(n^-2).  This 
quantity  HS^  is  then  assumed  to  be  an  estimate  of  the  value  of  the 
"average"  ir  within  laboratory  variance  ajjj  of  the  slope.  Another  quantity 
wh’ch  must  be  defined  is  the  between  laboratories  variance  of  the  slopes  B. 


Knj-lHBj-B)2 
b  ”  k-l 


In  (Reference  8)  B.  is  the  slope  obtained  for  data  of  i-th  laboratory 

Ai  * 

and  B  is  the  average  obtained  by  weighing  the  best  slope  for  each 
laboratory  as  follows: 


B  = 


X(nj-l)Bj 
?(n,— I) 


(9) 


It  is  assumed  that  the  best  composite  value  of  the  slope,  [B],  can  be 
determined  from  the  weighted  average  of  individual  slopes  defined  by: 
ZwjB| 

[B]  =  (10) 
i 

where  is  a  weighing  factor  which  is  further  defined  by  the  following 
relation: 

wj  =  [x  +  l/(n-n]-1 


where 


X  =  o-  5  /cr 


W 


(11) 


(12) 
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Tne  best  value  of  X  is  obtained  by  an  iterative  procedure  which  gives  a 
valu&of  G  =  0  for  the  following  relation: 


G  =Zt1(Br[B]}2-{k-i)o-2 


U3) 


The  first  estimate  of  X  and  thus  the  first  guess  of  the  weighing  factors, 

•5  2 

is  obtained  from  the  relations  given  for  a£  above  and  the  following 


relation  for  a 


2. 

fa¬ 


ff 


2 

b  = 


MSb~  MSyy 

"  ?M)2 

[^.(nj-D-  2{nr!)J/k‘l 

i 


(14) 


Subsequent  values  of  X  are  then  estimated  until  Equation  13  equals  0  to 
the  desired  approximation. 


The  best  value  of  (aHq)v ,  the  heat  vaporization  of  nickel  solid, 
which  was  obtained  by  utilizing  the  above  statistical  method,  was: 

(AHq)v  (s— v)=  102.9  ±  4.2kcol/mole. 

It  was  also  of  interest  to  obtain  the  value  of  B  for  a  given 
temperature  range  so  that  Equation  6  may  be  used  to  calculate  vapor 
pressures.  Since  the  data  utilized  covers  the  temperature  range  1273  to 
1709°K,  a  value  of  B  at  the  midpoint  of  this  range  was  calculated.  It 
should  be  noted  that  there  is  a  difference  in  specifying  the  average 
temperature  for  the  range,  1491°K,  as  the  midpoint,  and  the  average 
temperature  that  is  derived  from  the  reduction  of  a  set  of  data  by  the 
least  squares  technique.  The  value  of  B(1491°K)  (=AH^g-|/R)  is  given  by 
correcting  AHg  to  AH^g-j.  This  was  done  by  using  the  enthalpy  tables 
(Reference  16)  to  obtain  a  value  of  AH^g^  =  101.1  ±4.2  kcal/mole. 

There  is  a  maximum  error  of  0.65  kcal  over  the  temperature  range 
introduced  by  assuming  that  AH0  is  constant  over  the  range.  The  error 
of  0.65  kcal  is  the  systematic  error  at  the  end  points  of  the  range; 
to  a  good  approximation  the  total  error  over  the  range  may  be  estimated 
as  ±4.9  kcal.  Therefore,  the  average  value  AH°  (1273  to  1709°K)  is  given 
as  101.1  ±4.9  kcal/mole  and  the  value  of  B  =  50,884  ±2,466  over  the  range. 
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An  estimate  of  the  intercept  A  was  also  made  using  the  same  statistical 
treatment  which  was  used  in  deriving  the  value  of  B.  In  deriving  the 
value  A,  however,  the  data  of  Morris,  et  al  (Reference  11}  was  omitted 
since  this  data  refers  to  the  liquid  vaporization.  Furthermore,  since 
tne  parameters  in  Equation  6  must  relate  to  a  temperature  range  which 
does  not  involve  a  discontinuity  in  the  parameters,  such  as  caused  by  a 
phase  change  two  values  of  A  and  B  were  calculated  which  can  be  used  to 
calculate  the  vapor  pressure  over  the  whole  temperature  range  for  which 
data  is  available. 

The  "best  values"  from  the  second  law  determinations  for  A  and  B 
for  both  the  solid  and  liquid  ranges  are  given  by  Equations  15  to  18. 

For  the  reaction 

Ni{s)—  Ni(g)  AH§=  lOl.l  ±4.2kcal/mole  (15) 

and  the  vapor  pressure  over  the  temperature  range  1273  to  1709°K  is: 

In  Ps(Torr)  =  24.26  ±l.39-^8-8-^?^-6-  (16) 

For  reaction: 

Ni(  I)  —  Ni  (g)  Ah°8550  =  95.91  ±i  .06  kcal /mole  (17) 

there  is  only  one  value  of  AH0  which  was  obtained  from  data  in  Reference 
11.  This  data  was  statistically  treated  to  eliminate  pressure  data 
points  which  deviated  more  than  2a  (r.m.s.  deviation)  from  the  calculated 
values.  This  treatment  yielded  ths  following  values  of  A  =  22.222,  and 
B  =  48,271.  These  values  differ  slightly  from  the  values  of  A  =  22.314 
and  B  =  48,432,  reported  in  {Reference  11).  The  best  equation  for  the 
vapor  pressure  of  liquid  nickel  in  the  temperature  range  1816  to  1895°K 
is: 

In  Pf  (Torr)  =  22.222  ±0.29- £8i271_+533  (18) 
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In  subsequent  discussions  covering  the  third  law  "best  values," 
relations  between  A  and  the  free  energy  functions  will  be  discussed. 


TABLE  VI 

QUANTITIES  USED  IN  DETERMINING  SECOND  LAW  "BEST  VALUES" 
OF  AHq,  B  AND  A  FOR  VAPORIZATION  OF  NICKEL  (SOLID) 


Reference 

Number 
of  Data 
Points -n^ 

(“tyi 

kcal/mole 

Variance  (a) 

o\  (AHq) 

A 

Variance  (a) 
o\  (A) 

(c) 

24 

115.83 

5.48 

27.861 

3.389 

(4) 

18 

112.36 

1.56 

26.735 

0.144 

(9) 

10 

102.21 

0.34 

24.092 

0.0433 

(12) 

7 

93.99 

58.69 

21.354 

8.451 

(10) 

5 

89.36 

48.16 

20.955 

5.688 

(ID 

15 

103.10 

1.12 

(b) 

(b) 

(a)  Statistical  degrees  of  freedom  equal  (n--2)  for  each  determination. 

(b)  Liquid  nickel,  this  value  not  used  in  averaging. 

(c)  This  report  Part  I. 

(2)  Determination  of  Third  Law  Best  Values  of  B,  A  and  AHq 

Third  law  determinations  of  the  "best  valuec"  of  the  quantities 

B,  A  and  AHq  can  be  ■  ~de  by  utilizing  the  expressions  for  the  free  energy 

functions  and  the  values  of  these  functions  gi>'en  in  (Reference  3);  and 

the  data  reported  in  References  (4  through  12).  The  data  in  References 

5  through  7  were  omitted  because  it  was  inconsistent  with  other  data  as 

noted  above.  The  data  in  (Reference  10)  was  also  omitted  because  the 

value  of  AHq  obtained  from  this  set  of  data  deviated  by  more  than  2a 

from  the  average  AHq  obtained  by  simply  statistically  weighing  the  values 

of  each  set  of  data  by  the  number  of  data  points  reported  by  each 

laboratory.  Although  this  statistical  test  appears  to  be  somewhat  crude, 

it  should  be  remembered  that  such  an  averaging  process  overestimates  the 

variance  and,  therefore,  represents  a  weaker  test  of  nonhomogeneity  and 

2 

low  probability  of  the  validity  of  the  data  than  the  normally  used  x 
test  which  puts  an  upper  and  lower  bound  on  the  relative  values  of 
variances  for  rejecting  observations.  The  reason  that  only  an  upper 
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bound  was  used,  was  simply  due  to  the  fact  the  data  were  statist'cally 

2 

too  inhomogeneous  to  treat  by  the  x  test  which  utilizes  both  upper  and 
lower  bound. 


In  order  to  compare  the  equation  for  In  P  in  the  temperature  range 

1273  to  1709°K  with  that  obtained  by  means  of  the  second  law  as  noted 

above,  it  is  necessary  to  compare  equations  for  In  P  used  to  obtain 

third  law  values  of  AH0  with  that  used  to  obtain  the  second  law  values 

of  this  quantity.  In  the  use  of  the  second  law,  the  slope  is  equated 

to  the  quantity  (AH°/R)  where  f  is  the  temperature  corresponding  to  the 
T 

mean  inverse  temperature  as  noted  above.  Equation  5  used  to  obtain  the 
value  of  AHjjgg  by  the  third  law  which  can  be  expressed  as  follows: 

Rin  P  =  2  -  A-Hi.9-3  (19, 

g-s  \  T  }  T 

To  convert  Equation  19  to  the  "second  law"  form,  that  is  Equation  6 

In  P  =  A  +  •§•  (6) 

/AH£-AHSg8\ 

the  quantity  l — ^ - J  must  be  subtracted  from  and  added  to  the 

righthand  side  of  Equation  19  as  follows: 


Rln  p  _  1  _  ^G°~H298^  +  ^Ht~AH298j  _  AH298  ^^-^293^  (gp) 

This  equation  reduces  to  Equation  21.  The  tabulated  values  (References 
3  and  16)  of  (AHI-AHggg')  being  used  to  make  the  corrections 


Rln  P 


_  z_ 

g-s  \  T  / 


AHf 


(21) 
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( ST_H298\  *  ’s  a^n50St  constant  over 
As  will  be  shown  the  quantity  Z  -  ^ — j — j 

the  temperature  range  of  interest  so  that  its  average  can  be  taken  as 
the  constant  A  which  with  the  unit  of  pressure  in  Torr  is  given  by 

A  =  5=  /  ACpdT  +  In  760  <22) 

g-s  \  RT  /  RTJ  253 

where  ACp  =  Cp  (gas)  -  Cp  (solid),  and  is  a  function  of  T. 

The  slope,  B,  is  given  by  Equation  23 

a  =  -  J-P  ACpdT  (23) 

R  R-<298 

The  quantity  A(AHe)  is  a  constant  for  a  given  mean  temperature T. 


In  this  report  all  third  law  data  was  reduced  by  using  the  values  of 


given  in  Reference  3. 


An  estimate  of  the  error  in  the  tabulated  free  energy  functions  is 
limited  to  assuming  that  the  error  that  is  reported  by  Kelly  (Reference 
17)  for  enthalpies  in  this  temperature  range.  This  error  is  0.3%. 
However,  additional  errors  in  A  arise  due  to  tha  fact  that  Z 


-  ^GT~H2S8^ 


is  not  constant  over  the  temperature  range.  The  error  due  to  this 
variation  has  also  been  evaluated  by  finding  the  mean  value  of  this 
quantity  for  all  determinations  used  and  evaluating  the  r.m.s.  deviation; 
to  this  quantity  the  estimated  error  in  [o(AH°)]  was  added  to  obtain  the 
total  error. 


The  mean  value  of  A  for  the  temperature  range  of  interest  to  us  is 
that  at  mean  temperature  T  =  1535°K,  which  is  the  inverse  mean  temperature 
as  defined  on  p.  22.  This  method  of  obtaining  the  mean  temperature  is 
consistent  with  the  method  of  averaging  the  free  energy  functions, 
which  are  weighed  by  the  quantity  (1/T).  The  error  arising  in  the  use 
of  the  average  free  energy  function  over  temperature  range  1250°K  to 
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7710°K  is  not  random,  but  can  be  expressed  as  a  maximum  limit  of  the 
error  at  the  end  points  of  the  range.  The  maximum  error  over  the  range> 
including  the  0.32  error  in  the  free  energy  functions,  is  1.52  while 
the  average  deviation  is  0.72.  If  it  is  wrongly  assumed  that  the  errors 
are  random  and  the  r.m.s.  deviation  is  calculated,  then  an  estimated 
error  of  0.82  is  obtained.  The  above  is  a  statement  of  errors  which 
arise  in  the  use  of  the  average  value  of  S  =  23.50  at  T  =  1535°K.  The 
value  of  =  99.47  ±0.6  kcal/mole,  while  B  =  50,060  ±301.  The  best 

value  of  =  101.45  ±0.62  kcal/mole.  This  value  of  AHq  is  more 
reliable  than  the  second  law  value  of  102.9  ±4.2  kcal/mole  even  though 
they  are  nearly  equal.  The  near  equality  appears  to  be  fortuitous  when 
one  considers  the  range  of  the  data  used  in  the  two  determinations. 

The  best  third  law  equation  for  the  vapor  pressure  of  solid  nickel 
in  the  temperature  range  1250°K  to  1710°K  is  as  follows: 

In  Ps  (Torr)  =  23.50  ±0.35  -  -50’60°  ±30- 

The  third  law  parameters  of  equation  6  can  also  be  determined  for  liquid 
nickel  from  the  limited  amount  of  data  available  (Reference  11)  in  the 
temperature  range  1816  to  1895°K.  The  best  value  of  A  =  21.98  with  a 
mean  error  of  ±0.04  and  maximum  error  at  the  end  points  of  the 
temperature  range  of  ±0.05.  The  value  of  B  =  47,775  ±138  (mean  error) 
with  a  maximum  error  ±210.  As  pointed  out  above  there  are  systematic 
errors  in  assuming  that  A  and  B  are  constant  over  the  temperature  range 
so  that  the  maximum  errors  ’epresent  the  errors  at  the  end  points. 

The  third  law  equation  for  the  liquid  vapor  pressure  is  then  given  by: 

In  F/  ( Torr )  =  21 .98  ± 0.05  -  47»  77^-*2i-°- 
in  the  temperature  range  1816  to  1895°K. 
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An  interesting  aside  was  the  treatment  of  the  third  law  data  points 
to  obtain  a  composite  second  law  value  of  aHq  by  successively  eliminating 
data  points  whose  deviations  were  greater  than  2cf.  Starting  with  all  the 
original  points  reported,  i.e.  91  points, it  was  found  that  after  a  number 
of  trials  only  26  points  remained  which  were  from  the  following  sources: 
References  4,  8,  and  12.  A  value  of  aHq  =  99.3  ±0.38  Real /mole  was 
obtained.  This  calculation  further  indicated  the  lack  of  statistical 
homogeneity  in  the  composite  data. 


TABLE  VII 


THIRD  LAW  VALUES  OF  AH°  USED  TO  CALCULATE 


0 

"BEST  THIRD  LAW  VALUE" 

FOR  Ni(s)->Ni(g) 

o-(AH°) 

variance 

*H0 

kcal/mole 

1 v  0' 

Reference 

deviation 

kcal/mole 

o?(4H5) 

Part  I 

25 

102.63 

2.09 

4.368 

(4) 

16 

101 .29 

0.13 

0.0169 

(8) 

12 

100.74 

0.20 

0.04 

(12) 

7 

100.28 

0.84 

0.7056 

(11)  a 

16 

102.12 

0.07 

0.0049 

(5-7)a 

10 

73.60 

0.66 

0.4356 

(ioa) 

5 

96.76 

1.45 

2.103 

Best  Value 

101.45 

0.62 

0.38 

(a)  Not  used  in  determining  final  value  of  AHg. 

(b)  ni  -  number  of  determinations  used  after  eliminating  those  with 
deviation  >2a. 


Since  there  are  differences  between  the  statistical  treatments  of 
data  reported  herein  and  that  reported  by  Paule  and  Mandel  (References 
13  and  14)  a  discussion  of  these  differences  is  in  order.  The  differences 
in  the  statistical  treatments  arise  from  the  fact  the  data  available  on 
the  vapor  pressure  of  nickel  ^cks  statistical  homogeneity  and  thus  a 
composite  value  of  the  second  and  third  law  heats  of  vaporization  could 
not  be  generated  by  simply  treating  all  the  data  as  a  single  set  of  data 
and  obtaining  best  value.  The  treatment  which  was  used  take?,  into 
account  the  nonhomogeneity  of  the  data  by  statistically  weighing  the 
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results  of  the  individual  laboratories  which  were  used  in  generating  an 
over  all  "best  value".  The  statistical  test  used  for  rejecting  data  in 
this  report  is  a  deviation  of  the  observed  data  of  more  than  2o  (r.m.s. 
deviation).  Paule  and  Handel  (References  14  and  15)  used  as  F  test  which 
in  essence  put  limits  on  the  probability  between  0.975  and  0.025  on  the 
validity  of  the  data  they  used.  Although  a  2 a  error  limit  does  not  take 
into  account  the  number  of  degrees  of  freedom  for  each  curve  in  setting 
error  limits,  it  does  give  a  convenient  limit  to  the  maximum  deviation 
to  be  tolerated  which  has  been  used  by  other  authors  (References  21  and 
22). 


The  probability  for  the  value  of  t  (deviation/a)>2  was  between  5  and 
15%  with  a  large  proportion  of  the  cases  being  between  5  and  10% 

(Reference  22).  Therefore  in  general  this  validity  limit  is  not 
inconsistent  with  the  other  statistical  treatments  that  were  done.  It 
is  interesting  to  note  that  with  continued  elimination  of  data  points 
whose  deviations  are  greater  than  2c  there  is  a  reduction  in  the  number 
of  statistical  degrees  of  freedom  and  the  probability  that  a  given  value 
will  fall  outside  of  the  2a  limit  changes.  However,  it  is  also  important 
to  note  that  by  definition  the  homogeneity  (a  measure  of  self-consistency) 
of  a  given  data  point  is  measured  by  t  =  (deviation/a). 

SECTION  2  Thermodynamics  of  Nickel  Oxide  Vaporization 

In  the  course  of  reducing  the  available  data  on  the  vaporization  of 
NiO,  Grimley,  et  al  (Reference  4)  derived  a  value  for  heat  of 
dissociation  of  the  diatomic  molecules  NiO^j.  The  heat  of  dissociation 
thus  obtained  appears  to  add  consistency  to  the  values  of  heats  of 
vaporization  of  Ni^  which  were  reported.  Since  this  data  has  appeared, 
additional  data  has  been  reported  on  the  free  energies  and  enthalpies  of 
the  reactions  19-23,  which  form  a  closed  cycle. 
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References 

AHq  kcal/mole 

4 

Ni0(s)  ^  Ni0(g) 

124.4  ±0.4 

(19) 

4 

Ni°(g)—  Ni{g)+  0(g) 

92.7  ±0.9(oJ 

(20) 

(b) 

Ni(g)  -  Nits) 

—101.5  ±0.6 

(21) 

28 

°(gp  1/2  °2  tg) 

-58.969  ±  .005 

(22) 

23  to  27 

Nife)+  l/202~Ni0(s) 

-56.67^25)±.t3 

(23) 

(a)  Corrected  values  of  data  in  Reference  4.  Seep.  42of 
this  report. 

(b)  Best  composite  value.  See  p.  30. 

The  new  data  for  the  change  in  enthalpies  (reactions  19-23)  was  used 
in  conjunction  with  data  of  Reference  4  on  the  vapor  pressure  of  NiO^ 
to  obtain  the  best  value  of  the  heat  of  dissociation  of  NiO^,  The 
reported  data  (Reference  4'  oil  the  observed  partial  pressures  of  the 
various  gaseous  species  in  equilibrium  with  NiO^  were  corrected  by 
using  new  values  of  ionization  cross  sections,  Appendix  B,  and  new 
values  of  free  energy  functions  for  NiO^  and  NiO^j.  The  standard 
state  to  which  the  thermodynamic  functions  were  referred  was  0°K.  Since 
many  thermodynamic  functions  are  referred  to  298.13°K,  these  functions 
had  to  be  corrected  to  a  0°K  reference  state. 

A.  Thermodynamic  Functions  of  Solid  Nickel  Oxide  (NiO)$ 

The  va1_~  of  (H°298  13  "  hq)  for  Ni0(s)  was  obtained  using  the  heat 
capacity  data  reported  by  King  (Reference  29).  This  data  is  given  in 
Table  III.  A  calculation  was  made  to  determine  the  difference  in 
enthalpy  of  (NiO) s  between  54.28°K  and  0°K.  This  was  done  by  extrapolating 
the  data  in  Reference  4  by  the  use  of  the  following  expression  for  the 
enthalpies: 


H|  -  Hq  =d(^)  +  e(^P),  where  o(-^p)  (References  30  and  31)  is  the 

Debye-Enthalpy  Function  and  E  (^-)  »  (Reference  31),  the  Einstein- 
Enthalpy  Function.  The  number  (404)  is  an  apparent  value  of  the  Debye 

33 


l  1 1  -.I  r  •jfttiBhi'irtairi&iV'i.i.iii  ■ 


AFAL-TR-72-217 


temperature  and  (620)  is  an  apparent  value  of  the  Einstein  temperature. 

This  calculation  gave  a  value  for  (H°5428  '  H(P  =  11,56  ^Z^16, 

The  quantity  (H°g8  13  -  H°4 2Q )  was  calculaced  by  adding  increments  of 

enthalpy  over  limited  temperature  intervals.  Taking  the  average  value 
of  the  heat  capacity,  C  ,  over  the  interval  defined  by  points  (n+1)  and 
(n),  where  n  represented  a  point  in  the  numerical  sequence,  the  enthalpy 
difference  between  T  and  Tfi  was: 

-  »Tn)  .Krtt  T„  )  (24) 

The  sum  of  the  enthalpy  increments  between  any  two  temperatures  gives 
the  enthalpy  difference  between  the  temperatures.  This  procedure  as  noted 
gave(H°g8  -  Hg)  =  1609.58  cal/mole.  i.e. 

(u°  -  \\°A  oq)  =  1598.02  cal/mole;  adding  to  this  the 

quantity  (H°54  28  -  H°)  =  11.56  cal/mole,  a  value  ot  (H298J3  -  H0)  - 
1609.58  cal /mole  was  obtained. 

The  free  energy  functions  of  Ni°(s) :  ["(gt"H0^T3  i,e* 

-(G°-H°)/T  =  Sj-(Hj-Hq)/T, 

were  calculated  from  tabulated  and  extrapolated  values  of  standard 
entropy  Sy,  at  the  temperature  T,  and  the  enthalpy  difference  expressed 
as  (H°-H°)/T  (Reference  18).  Since  enthalpy  data  is  tabulated  as  (Hy-H  2gg 
correction  as  detailed  above  was  applied  giving  the  thermodynamic  data 
used,  as  shown  in  Table  IX. 

B.  Calculation  of  Thermodynamic  Functions  of  NiO^ 

The  thermodynamic  functions  of  (NiO)g  were  calculated  by  assuming 
that  the  NiO  gas  molecule  is  a  perfect  gas  composed  of  diatomic  rigid- 
rotor-oscillators  as  in  the  case  of  Ni'2  (Appendix  A).  The  contributions 
of  the  individual  degrees  of  freedom  of  the  possible  motions  of  the 
molecule  to  the  free  energy  function  were  determined.  The  contribution 
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TABLE  VIII 


T,  °K 

c 

S 

cal/°K/csole 

HEAT  CAPACITIES  OF  Hi0(sj 

C„,  cal/°K/nole 

AT  "U  P 

(4H3  -Hf  ) 
‘n+1  'n 
cal/raole 

54.28 

0.900 

4.65 

1.006 

4.678 

58.93 

1.112 

4.56 

1.227 

5.595 

63.49 

1 .341 

4.61 

1.462 

6.740 

68.10 

1.582 

4.37 

1.700 

7.429 

72.47 

1.818 

4.27 

1.935 

8.262 

76.74 

2.051 

3.48 

2.150 

7.482 

80.22 

2.248 

3.59 

2.346 

8.422 

83.81 

2.444 

10.85 

2.750 

29.84 

94.66 

3.056 

10.31 

3.348 

34.52 

104.97 

3.639 

9.60 

3.899 

37.43 

114.57 

4.158 

10.03 

4.425 

44.38 

124.60 

4.691 

11.28 

4.980 

56.17 

135.88 

5.269 

9.67 

5.503 

53.21 

145.55 

5.737 

10.31 

5.979 

61.64 

155.86 

6.221 

9.92 

6.430 

63.79 

165.78 

6.639 

10.06 

6.840 

68.81 

175.84 

7.041 

9.94 

7.234 

71.91 

185.78 

7.427 

10.16 

7.606 

77.28 

195.94 

7.785 

10.19 

7.960 

81.11 

206.13 

8.134 

10.06 

8.296 

83.46 

216.19 

8.457 

9.66 

8.610 

83.17 

225.85 

£.762 

10.22 

8.906 

91.02 

236.07 

9.050 

9.56 

9.178 

87.74 

245.63 

9.306 

10.66 

9.452 

100.76 

256.29 

9.597 

9.86 

9.721 

95.85 

266.15 

9.845 

9.86 

9.963 

98.24 

276.01 

10.08 

10.42 

10.20 

106.28 

286.43 

10.32 

9.51 

10.44 

99.28 

295.94 

10.55 

2.21 

10.55 

23.52 

H°298.rH°54.28  = 

1598.02 

LJO  MO  _ 

H54.3'H0 

11.56(b) 

H298.13"H0 

1609.58  cal /mole 

(a)  Reference  29. 

(b)  Calculated  by  using  the  following  expression  for  the  enthalpy 
utilizing  Debye  functions  and  Einstein  functions: 


7422) 

,  f/585\ 

\54.3/ 

(c)  Johnston  and  Marshall  (Reference  8)  report  H°QQ  =  1630  cal, 
using  data  H.  Seitz, 
expression  for 


J.  OeWitt,  H.J.  McDona1d  (Reference  33)  with 
,  between  0°K  and  68°K. 
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FREE  ENERGY  FUNCTIONS  FOR  NiO(sj  (REFERENCE  18) 


-(Gf-H5/T)s 

T°K 

cal/°K/mole 

1500 

18.43 

1600 

19.19 

1700 

19.91 

1800 

20.61 

to  the  free  energy  functions  by  the  translational  degrees  of  freedom, 
~[Gf-H§/T]t,  is  given  by  Equation  41.  The  rotational  partition  function 
was  determined  by  utilizing  Equation  43,  and  noting  that  the  symmetry 
number  a  -  1,  and  that  the  rotational  constant  has  the  value  given  in 
Table  IX.  The  electronic  contribution  is  given  by  Equation  42.  The 
electronic  partition  function  required  for  Equation  42  was  determined  as 
noted  in  the  next  paragraph  utilizing  the  energy  levels  given  in  Table  X. 

(1)  Estimate  of  Electronic  Contribution  to  (NiO)g  Partion  Function 

In  determining  the  electronic  contribution  to  the  partition 
function  of  NiO  molecule,  and  subsequently  to  the  thermodynamic  functions, 
the  degeneracy  (statistical  weight)  of  the  spectroscopic  state  (electronic 
energy  level)  must  be  determined.  Previous  estimates  of  the  partition 
functions  were  incorrect,  because  the  spectroscopic  states  of  NiO  were 
considered  to  be  derived  from  the  nickelous  ion  (Ni++)  (Reference  36), 
therefore  it  is  necessary  to  redescribe  the  atomic  and  molecular  binding 
in  the  diatomic  molecule  NiO  in  order  to  re-estimate  the  partition 
function. 

The  assumption  made  ir  Retv.-.-ero>>  36  that  the  (NiO)  molecule  is  "'onic 
ignored  the  preponderance  of  the  evidence  on  the  electronic  characteristics 
of  diatomic  molecules  which  indicates  that  most  diatomic  molecules  do 
not  exist  as  ionic  molecules  (Reference  37).  Therefore,  the  description 
of  the  electronic  states  of  the  NiO  diatomic  molecule  which  must  be  used 
ere  those  consistent  with  the  eU<-tronic  states  of  atomic  molecules. 

To  test  the  idea  that  (NiO)  is  atomic,  the  instability  of  the  -onic 
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1 .40  ev  (Reference  40).  The  molecule  [Ni  0  ]  Is  even  more  unstable 
then  [Ni+0~]  since  I(Ni+)  =  18.2  ev  (Reference  28)  and  E(0“)  =  zero 
(Referemce  35). 

Thus  the  total  difference  in  the  quantity  [I-E]  in  gei’g  to  the 
second  ionization  potential  of  Ni,  to  form  [Ni‘  ]  and  to  :he  second 
electron  affinity  of  oxygen  to  form  the  [0~] ,  species  is  [25.81  -  1.465] 
ev  or  24.34  ev.  This  difference  is  obviously  too  large  for  the  molecule 
[Ni.  0-]  to  be  stable  with  respect  co  the  covalent  NiO(g)  molecule. 

To  test  further  the  validity  of  the  atomic  molecular  structure  for 
[Ni 0]  the  ionicity  of  the  [NiO]  molecule  was  estimated  by  taking  the 
difference  between  the  average  heats  of  dissociation  of  the  diatomic 
molecules  of  the  elements  and  the  heat  of  dissociation  of  gaseous  NiO. 

The  heat  dissociation  used  for  Nig  was  that  found  by  recalculating 
Kant's  data,  (Reference  42)  i.e.  Dg(Nig)  =  61  kcal/mole  (p.  52).  The 
heat  of  dissociation  of  0g  used  was  that  used  whenever  electronegativities 
are  used  in  conjunction  with  thermodynamic  data  and  is  equal  to  the  heat 
of  dissociation  of  two  single  bonds,  i.e.  66.5  kcal  instead  of  the  measured 
Dg  =  118  kcal/mole  (Reference  45).  The  heat  of  dissociation  of  NiO  was 
estimated  as  85  kcal  (Reference  4)-  Defining  A,  as  a  measure  of  the 
ionicity,  one  finds  that  for  [NiO]: 

A  =  1/2  D^tC^)  +  1/2  Ni2)  -  Dq(NiO)  =  21k  col, 
or  that  the  ionicity  is  relatively  low. 

Having  concluded  that  ionicity  of  NiO  is  low  one  is  lead  to 
assumptions  concerning  the  nature  of  the  ground  electronic  state  of 
[NiO].  Such  assumptions  can  be  made  on  the  basis  of  the  observed 
electronic  transitions.  However,  since  these  transitions  have  not  been 
classified  for  NiO,  one  must  guess  at  the  electronic  configuration  of 
the  N"'0  molecule.  The  ground  state  could  be  classified  either  as  a  ^ 
or  ^A  configuration  since  these  are  possible  electronic  configurations 
of  combined  [0]  and  [Ni]  atoms,  which  have  been  determined  by  taking  the 
possible  combination  of  atomic  states  as  given  by  Herzberg  (Reference  44). 
The  jssumption  that  the  ground  state  should  be  a  state  leads  to  the 
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assumption  that  the  excited  states  should  he  designated  as  states. 

This  labeling  is  the  most  consistent  with  our  ignorance  and  information, 
and  it  does  not  appear  to  violate  any  selection  rules  nor  is  it  contrary 
to  a:*y  spectroscopic  observations.  Having  re-estimated  the  values  of  the 
electronic  partition  functions  by  using  the  expressions  given  in  Table 
XIII;  other  available  molecular  data  was  used  to  recalculate  the  free 
energy  function  of  NiO  as  illustrated  in  Table  X. 

TABLE  XI 


CORRECTED  PRESSURES  OF  Ni0^gj  AND 
THIRD  LAW  HEATS  OF  VAPORIZATION:  NiO$+NiOg(a) 


Pressure^ 

AH  § 

T°K 

atm 

kcal/mole 

1575 

28.59x10"]° 

56.92x10"^ 

53.46x10"q 

125.35 

1587 

124.09 

1596 

124.97 

1606 

8.52x10~q 

124.22 

1624 

1C.90x10~q 

124.76 

1625 

11. 35x1 0"X 

124.*’0 

1630 

13.91x10~q 

124. 10 

1646 

20. 70x1 Oq 

125.16 

1651 

23.53x10  q 

124.29 

1657 

24.68x10  q 

124.32 

1659 

28.52x10  q 

124.14 

1673 

37.82x10  q 

124.20 

1673 

36.22x10 "X 

124.35 

1679 

42.18x10”q 

124.30 

1679 

46.86x10 

123.93 

1684 

45.26x10  q 

124.38 

1684 

47 . 62xl0~o 

7.69x10  q 

124.21 

1707 

124.20/.  x 
125.68'°' 

1709 

51.92xl0'y 

Pressures  =  0.641 

x  Pressures  reported  in 

Reference  4. 

Omitted  from  final 

determination  of  AHg  as  shown  on  p.  41. 
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C.  Examination  of  the  Data  on  the  Enthalpy  of  Vaporization  of  NiO^j 

Since  the  only  available  data  on  the  partial  pressures  of  NiO^j  in 
equilibrium  with  ^0(s)  is  that  reported  in  Reference  4,  this  data  was 
examined,  and  as  noted  in  Appendix  B-l,  the  reported  pressures  were 
corrected  for  the  latest  data  available  on  ionization  cross  sections. 

The  corrected  data  on  the  equilibrium  partial  pressure  of  NiG^j  tor 
Equation  19  is  given  in  Table  XI  as  well  as  the  heats  of  vaporization 
calculated  by  the  third  law.  Utilizing  the  statistical  criteria  as  given 
in  Part  II,  Sec  B1  and  B2,  a  best  value  of  AHg  =  124.4  ±0.4  kcal  was 
obtained.  By  utilizing  the  relations  indicated  in  Part  II,  Sec  B2, 
values  of  the  parameters  A  and  B  for  the  vapor  pressure  equations  were 
also  calculated  by  averaging  the  values  of  the  free  energy  functions, 
over  the  temperature  range,  1575  to  1707°K,  taking  a  value  of  T  =  1646°K, 
and  correcting  the  third  law  values  of  AHg  to  AH^g.  The  equation  for 
the  vapor  pressures  of  NiO^,  in  equilibrium  with  NiO^,  can  thus  be 
expressed  as 

ln  pTorr(Ni0)g=  25  52  ±a05  "  (25) 

The  error  limits  were  set  by  the  observed  error  in  AHg  plus  the  reported 
error  cf  0.45£  in  the  enthalpy  of  solid  (Reference  17).  The  error  in  A 
is  systematic  since  it  is  temperature  dependent;  at  the  extremes  of  the 
temperature  range,  1575°K,  it  is  equal  to  +0.16,  and  at  1709°K  it  is 
equal  to  -0.13. 

A  redetermination  of  the  second  law,  AHg,  of  reaction  (Reference  19) 
was  also  performed.  The  best  least  squares  value  obtained  for  AH^g  = 
127.6  ±3.1  kcal/mole.  The  best  value  of  AKg  =  133.3  ±3.2  kcal  per  mole. 
The  best  value  of  the  parameters  A  and  B  are  shown  in  the  vapor  pressure 
equation  (Reference  26). 


In  pTofr(NiCV  2792  ±0-96“ 


(26) 
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The  discrepancies  between  the  second  and  third  law  values  of  the 
parameters  in  the  vapor  pressure  equation  appear  similar  to  those  obtained 
for  nickel.  Part  II,  Sec  B1  and  B2.  Due  to  the  greater  consistency  of 
the  third  law  determinations,  the  third  law  values  were  adopted  as  the 
best  values.  Tables  X  and  XI  show  the  data  used  in  determining  the  tfird 
law  parameters.  Gri ml ey (Reference  4)  reported  value  of  AHq  =  129.5 
±5  kcal/mole.  The  vapor  pressure  data  of  Johnston  and  Marshall  (Reference  8) 
on  NiO,  was  not  utilized  in  this  ca^e  because  there  was  diff>'  ilty  in 
differentiating  between  the  loss  of  weight  from  the  sample  due  to  NiO^j 
and  that  due  to  Ni^,  leaving  doubt  as  to  the  accuracy  of  the  results. 

D.  Examination  of  the  Data  on  the  Heat  of  Formation  of  NiO^ 

To  complete  the  evaluation  of  the  available  information  of  the 
enthalpy  changes  for  Reactions  19  to  23  an  examination  of  the  available 
data  on  the  heats  of  formation  NiO^  (Reaction  23)  was  performed.  In 
the  process  of  performing  this  evaluation,  the  mass  spectrometric  data  of 
Grimley,  et  al  (Reference  4)  on  the  partial  pressure  of  02  was  corrected  j 

by  utilizing  new  information  on  the  ionization  cross  sections  of  02  and  j 

0  (see  Appendix  B-.l).  ! 

Table  XII,  lists  the  0?  pressures  reported  by  Grimley  corrected  for  ( 

errors  in  the  ionization  cross  sections  as  noted  in  p.  67.  On  the  basis 
of  this  data  and  utilizing  the  second  and  third  law  methods  for  determining 
heats  of  dissociation  of  NiO^  as  per  Equation  23,  the  values  obtained 
were  as  follows:  AHq  =  57.37  ±0.30  kcal/mole;  AH2g8  =  57.94  ±0.30  kcal/ 
mole.  The  best  literature  value  of  AH2g8  -  57.240  ±.13  kcal;  and  that  of 
AHq  =  56.67  ±0.13  kcal/mole.  Reference  26. 

It  is  of  interest  to  note  that  the  r.m.s  deviation  of  Grimley's 
corrected  values  of  the  oxygen  pressures  from  the  best  values  of  these 
pressures  reported  in  Reference  26  is  42%.  The  comparison  was  made  by 
interpolating  the  pressures  reported  in  Reference  26  for  the  temperatures 
reported  by  Grimley  (Reference  4).  This  relatively  good  agreement  between 
the  mass  spectrometric  technique  and  the  high  temperature  E.M.F.  techniques 
utilized  in  obtaining  the  data  analyzed  by  Kellogg (Reference  26),gives 
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TABLE  XII 

CORRECTED  DATA  FOR  REFERENCE  4  02  PRESSURES 
AND  THIRD  LAW  AH^FOR  REACTION:  Ni^  +  1/2  02  -»■  NiO^ 


Temperature 

Pressure  oJa^ 

AH§ 

°K 

Torr 

d 

kcal/mole 

1575 

1.74  x 

10"4 

-57.87 

1587 

2.58 

II 

-57.65 

1596 

2.76 

II 

-57.85 

1606 

3.31 

II 

-57.93 

1624 

7.43 

II 

-57.23 

1625 

6.22 

II 

-57.55 

1630 

6.98 

II 

-57.52 

1646 

1.02  x 

10“3 

-57.47 

1651 

1.24 

II 

-57.24 

1657 

1.64 

II 

-57.04 

1659 

1.44 

II 

-57.29 

1673 

1.88 

II 

-57.31 

1673 

1.95 

II 

-57.26 

1679 

2.12 

II 

-57.34 

1679 

2.47 

II 

-57.08 

1684 

2.58 

It 

-57.16 

1684 

3.14 

II 

-56.83 

1707 

4.01 

II 

-57.17 

1709 

4.02 

It 

-57.25 

(AH°)f  (NiO)  =-57.37  +0.30  kcal/mole. 

(a)  Pressure  of  0?  =  1.25  Pressure  0?  in  Reference  4. 
see  p.67. 


confidence  that  the  mass  spectrometer  can  be  used  to  determine  quantitatively 
the  flux  from  a  Knudsen  cell  and  thus  the  pressure  inside  the  cell. 

Determination  of  Og(NiOg) 

The  determination  of  D° ( N i 0 ) g  as  per  Equation  20  is  now  reduced  to 
completing  the  thermodynamic  cycle  defined  by  Equations  19  through  23. 

The  best  value  of  Dg(NiO)  =  92.7  ±0.9  kcal/mole.  The  error  in  Dg  being 
determined  as  the  r.m.s.  of  the  sum  of  the  squares  of  percentage  errors 
of  the  quantities  in  the  cycle. 
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It  should  be  noted  that  although  the  value  of  Dq  (HiO)  =  S2.7  kcal/mole 
differs  from  the  value  of  Dq  =  S6.5  kcal  reported  in  Reference  4,  it  is 
still  consistent  with  this  and  other  reported  values:  9q  =  87  kcal 
suggested  by  Trivedi  (Reference  45)  and  Dq  <4.2  ev  (96.6  kcal/mole) 
suggested  by  Kuldt  and  Lagerquist  (Reference  35). 
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PART  III 
APPENDIX  I 

DETERMINATION  OF  D°  FOR  Nig 

It  was  necessary  to  obtain  a  value  of  the  dissociation  energy,  Dg, 
for  the  Ni0  molecule  in  order  to  justify  the  assumption  that  the  ground 
state  of  NiO  is  a  Z  state.  A  survey  of  the  literature  indicated  that 
Kant (Reference  42)  had  made  the  only  reported  determination  of  Dj^Nig). 

The  data  reported  by  Kant  was  critically  examined,  and  several  errors 
were  noted  in  the  data  reduction,  therefore,  the  reported  data  had  to  be 
reconciled  with  the  reported  results.  For  example,  Kant (Reference  42) 
reported  that  the  experimentally  determined  ratios  of  the  mass  numbers 
which  he  stated  contributed  most  significantly  to  the  composition  of  the 
gas  phase  Nig^Ni^16,  Nig18,  and  Ni^28)  were  as  follows  116:118:120  =  1.0: 
0.77:0.15,  and  he  further  stated  that  these  ratios  corresponded  to  that 
expected  from  known  isotopic  distribution.  As  shown  in  the  following 
analysis  these  ratios  are  not  what  is  expected  from  isotopic  distribution 
in  the  natural  abundance.  An  examination  of  the  natural  isotopic 

to 

abundance  of  nickel  gives  the  following  isotopic  distribution:  Ni  - 
67.88%;  Ni60  -  26.32%;  Ni61  -  1.19%;  Ni52  -  3.66%;  and  Ni63  -  1.08% 
(Reference  46).  These  ratios  in  turn  can  be  used  to  determine  the  ratio  of 
diatomic  molecules  of  various  isotopic  composition  either  by  statistical 
determination  of  the  distribution,  or  from  the  ratio  of  the  equilibrium 
constants  for  the  reactions  for  the  dissociation  of  the  diatomic  molecules 
such  as: 

Ni'  +  Nil  Ni|i  (27) 

The  equilibrium  constants  for  these  reactions  differ  from  each  other 
by  the  symmetry  numbers  where  =  2  for  homonuclear  molecules, 
and  =  1  for  heteronuclear  molecules.  The  equilibrium  constant  for 
the  reactions  involving  the  isotopes  i,j  are  equal  to  (K/o - • )  where  K  is 

•  J 

the  equilibrium  constant  for  reactions  involving  isotopes  for  which  i  £  j 
(References  47  and  48).  Equilibrium  constants  and  reactions  for  all  the 
isotopic  combinations  of  i  and  j  can  be  written  by  simply  rotating  i  and 
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j,  over  all  the  possible  values  of  i  and  ,j.  In  applying  the  relations 

for  the  equilibrium  constants  as  given  in  (Reference  28),  it  is  noted 

that  free  atom  ratios  Ni^Ni^iNi^’.Ni^Ni63  are  the  same  as  in  the 

natural  isotopic  abundance  of  nickel;  therefore,  the  ratios  such  as 
i  i  i  k 

[Ni^/Nig  ]  can  be  determined  from  the  ratios  of  the  equilibrium  constants 
for  the  individual  isotopic  reactions. 


K  _  [n$] 

<r,I  "  [Nii][Nii] 

°~ik  _  (  Ni2* ]  (  Nii] 
a\i  ‘  'Mjk' 'Ni1' 'Nil ' 

fjk  _  /Ni2*WNiH\ 
a\\  'Ni^k/  V  Nii  / 


(28) 

(29) 

(30) 


For  example:  in  the  above  equations,  letting  i  =  60,  and  k  =  58  then  if 
i  =  j  it  follows  that  i,  j  =  i»  i  =  60,  60,  while  i,  k  =  60,  58; 

CO  CQ 

furthermore,  =  2  and  =  1.  respectively,  and  the  ratio  Hi  :Ni  = 
67.52:26.32  =  2.565  and  one  obtains  the  following  value  for  the  ratio  of 
the  diatomic  molecules. 


[Nif°>60/Ni|8«60]  =  1/2  [Ni60/Ni58]  =  |.283 


Since  the  mass  spectrometer  determines  mass  numbers,  all  possible 
combinations  of  i  and  j  whose  sums  are  equal  to  a  given  mass  number 
contribute  to  the  current  for  a  given  mass  number  and,  therefore,  must 
be  included  in  determining  the  relative  ratios  of  the  mass  numbers.  For 

example  and  Nil^3  can  be  formed  only  by  the  combination  Ni33Ni^3 

59  59  .  120 

and  Ni  Ni  ,  respectively,  while  the  molecule  Ni 2  can  be  formed  from 

the  combination  Ni®°Ni6®  and  Ni53Ni62. 


Computation  of  the  ratios  of  all  possible  combinations  of  isotopes 
which  contribute  to  abundance  Ni^  ,  Ni^  >  and  Ni^  *  yield  the  following 
ratios:  [Ni^iNi^^Ni^]  =  1:0.78:0.257.  The  discrepancy  of  these 
ratios  to  those  reported  by  Kant (Reference  42),  given  above,  puts  the 
experimental  accuracy  of  Kant's  data  under  suspicion  but  because  the  data 
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presented  Is  independent  of  his  determination  of  the  molecular  isotopic 
ratios,  the  dissociation  energy  [D^Kig)]  can  still  be  estimated  by 
using  Kant's  data,  and  the  relations  he  used.  The  relation  for  Dg 
(Equation  32)  which  was  derived  by  Drowart  and  Honig  (Reference  49), 
appears  to  have  an  error  in  the  constant,  therefore,  a  derivation  will 
be  given  of  Equation  32  so  that  the  value  of  the  constant  can  be  checked 
by  the  reader. 

Relations  for  Determining  Dg(Ni2)  from  Spectroscopic  Data 

Equation  32  is  the  relation  given  by  Drowart  and  Honig (Reference  49) 
for  the  dissociation  energy  of  a  diatomic  molecule  in  terms  of 
spectroscopically  derived  quantities.  In  this  particular  case  it  is 
assumed  that  the  diatomic  molecule  is  the  dimer  of  the  atom. 

2~30^KT  =  ~  log  pl  *mfI^  +  log  “ log  ^y2/y|) 

4-  3/2  logT  ir  log  Mj  +  log  Q^e| 

—  log  Q2el  “  2  logr2“  log  Q2v  +  3.2771  (32) 

In  Equation  32,  P^Torr)  is  the  pressure  of  the  monomer  in  Torr, 
(Ig/I-j).  is  the  ratio  of  diatomic  current  to  atomic  currents.  The 
diatomic  current  is  of  course  the  sum  of  the  currents  of  the  isotopically 
different  masses  which  comprise  the  dimer  and  the  atomic  current  is  the 
sum  of  isotopic  atomic  currents.  It  is  further  assumed  that  (I^/lj)  = 
[P2/P^]  where  P2  is  the  pressure  of  the  diatomic  molecule.  The  symbols 
y-|  and  Y2  represent  the  products  of  ionization  cross  section  and  multiplier 
efficiency  for  the  atom  ion  and  diatomic  ion-molecule  respectively; 

T  is  the  temperature  in  °K;  M-j  is  the  atomic  weight  of  the  atom,  which 
in  the  case  of  Ni  equals  58.7.  Q2e^  and  are  electronic  partition 
functions  for  Ni2  and  Ni,  respectively.  Q2y  is  the  vibrational  partition 
function  for  (Ni2);  r2  is  the  radius  of  the  diatomic  molecule,  (Ni2); 

R  is  the  gas  constant,  and  3.2771  is  a  constant.  The  quantities  used 
by  Kant  (Reference  49)  and  those  used  herein  to  recalculate  the  values 
of  Dq(N1 2)  are  given  in  Table  XIII. 
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TABLE  XIII 

to 


DATA  USED  TO  CALCULATE  Dg  BY  EQUATION  32 


Term 


CY2(Ni2)/^Ni)3 


H, 


Qel(Ni2) 


Expression  or  Value 

1  (a) 

58.71 ^  g/mole 

Ab) 


qe1(Ni)c-d 

9  +  7  exp 

[.MZ; 

j  +  5  exp 

m 

+  7  exp 

J  2948 j 

+  5  exp 

m 

+  3  exp 

[_  2466] 

+  5  exp 

r  49091 
■  T  J 

r 

2.30A^ 

n  (a),  (e) 

Vv 

[l  -  exp 

467. 7J 

-1 

constant  in  relation  32 

3.2771;  2.8247(a) 

(a)  Value  given  in  Reference  42. 

(b)  This  value  of  Q£] (Nig)  =  1  is  obtained  by  assuming  that  the 

ground  state  of  (Nio)  is  ^E.  This  differs  from  the  assumption 

■? 

in  Reference  42.  The  choice  of  the  E  g.s.  is  based  on  the 
possible  electronic  configurations  of  Nig  [Reference  37]. 

(c)  Reference  42  has  a  typographical  error  in  this  expression. 

(d)  The  electronic  levels  (cm-^),  and  their  degeneracies,  of  the 
nickel  atom  which  contribute  to  the  partition  function  are  as 
follows:  0.0(9),  204.8(7),  879.8(5),  1332.2(7),  1713.1(3), 
2216.5(5),  3409.9(5). 

(e)  The  vibrational  frequency  of  Nig  was  assumed  to  be  325  cm"^. 

The  value  was  obtained  by  taking  1.25  x  Debye  frequency  of  the 
solid  as  in  Reference  42. 
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The  equations  and  numerical  constants  used  in  deriving  Equation  32 
are  obtained  from  the  relations  given  by  Herzberg  (Reference  50)  and  the 
OANAF  Tables  (Reference  28) .  The  free  energy  functions  which  are  utilized 
are  those  for  the  species  in  Reactions  33,  34,  and  35. 


Ni2(g)-2Ni(g) 

ahq  =  d§ 

(33) 

2Ni(s)-Ni2(g) 

(AH^)v(Ni2) 

(34) 

Ni{s)  —  Ni(g) 

(AHo)v(Ni) 

(35) 

The  changes  in  free  energy  functions  for  these  reactions  are  related  to 
equilibrium  constants  for  the  reactions  and  may  be  combined  as  shown  in 
Equation  36  and  37  for  Reaction  33. 


In  Equation  37  the  relations  between  the  change  in  standard  free  energy 
AG°,  and  the  equilibrium  constant  was  utilized;  as  well  as  the  expression 
for  the  equilibrium  constant  in  terms  of  the  partial  pressures  and 
activities  of  the  species  involved. 

Rearranging  terms  and  combining  Equations  36  and  37  the  following 
equations  are  obtained: 

109  [kp]„  33  +  {§  =  «8) 

,„g  Rim)  + 1..,  +  *  *  - +(^$M 

Finally 

-0—  =  -log  P  (Ni)  +  log  - 2  (£^dfi)[r  -(g)]  [Nijg)l  (40) 

2.303RT  P(Ni)  '  RT  /L  J  \  RT  / L  Z  J 
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The  quantity 


S°-Hq 


[Ni(g)3  is  the  free  energy  function  of  nickel 


atom  and  can  be  calculated  from  the  sum  of  the  relations  for  the  electronic 
and  translational  contributions  to  the  free  energy. 


The  translational  contribution  to  the  free  energy  function  is  given 
by 

-(ffrgo)  trons  =  6.863753  logM  +  11.439588  log T-  7. 283739  col/deg/mole  (41) 

In  the  above  equation,  M  is  the  molecular  weight  of  the  species  and  T  is 
the  temperature  in  °K. 


The  electronic  contribution  is  given  by: 

el  =4.575835  log?g-,exp-  cgi/deg/p^e  (42) 

where  is  the  energy  of  the  ith  electronic  level,  g.-  is  the  degeneracy 
of  ith  level  and  summation,  2^  is  taken  over  all  electronic  levels. 

For  diatomic  molecules  the  translational  and  electronic  contribution 
to  the  free  energy  function  is  given  by  Equations  41  and  42  and  vibrational 
and  rotational  contributions  are  given  as  follows: 


-^•yS.jrot  =  -4.575835  log  +  0.953116  (B/T)  +  0.0457127  (B/T)2 

-0.723040  (43) 

In  the  above  expression  o  is  the  symmetry  number  which  is  equal  to 
2  for  homonuclear  molecules  and  1  for  heteronuclear  molecules.  When 
considering  a  mixture  of  isotopes,  the  symmetry  number  for  diatomic 
molecules  of  the  same  element  is  equal  to  2.  However,  when  the 
distribution  of  isotopic  mixtures  of  diatomic  molecules  of  the  same 
element  is  being  determined  the  molecules  have  a  symmetry  of  1  when  the 
molecules  are  composed  of  different  isotopes  and  2  when  they  are  composed 
of  the  same  isotope  (Reference  47). 
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The  rotational  constant  B  is  given  by  the  relation 

B  _  2.798890  x  I  O'59  {44) 

I±t2 

where  r  is  the  radius  of  the  diatomic  molecule  *  y  is  the  ^educed  mass 
which  for  homonuclear  molecules  is  equal  to  (M^/2).  Where  is  the 
atomic  weight  of  the  element.  When  the  value  of  B  is  obtained  from 
spectroscopic  data  it  is  given  by  B  =  (Bg  where  B  is  the  rotational 
constant  measured  for  a  given  vibrational  state,  Bft  is  the  value  of  B 
at  the  equilibrium  distance  between  the  atoms,  and  ctg  is  a  correction  to 
be  applied  due  to  change  of  the  moment  of  inertia  with  vibrational  energ-'. 

If  it  is  assumed,  as  it  is  in  most  cases,  that  the  rotator  is  rigid,  then 
the  value  of  B  to  be  used  is  Bg.  In  principle  however,  B  must  be 
determined  for  each  vibrational  state,  but  since  the  variation  of  B  with 
vibrational  energy  is  small  even  for  high  temperatures,  the  rigid  rotator 
vibrator  approximation  is  accurate  enough  for  most  thermodynamic  applications. 
The  vibrational  contribution  to  the  free  energy  function  is  given  by: 

_(6°-HQ)uih  -  -4.575835  log  [  I  -exp-^j  col  /  deg /mole  (45) 

where  y  =  (1.43840  w/T)  and  w  is  the  vibrational  frequency  in  cm  . 


Summing  the  expressions  for  each  degree  of  freedom  for  the  Individual 

f  G°-H‘  1 

species  one  obtains  the  value  of  A  [ — T-^  J  ^or  dissociation  reactii 
of  a  homonuclear  diatomic  molecule  such  as  Nig.  For  example,. evaluating 
some  of  the  terms  such  as: 


§£  .  2  x  2.798890  x  IQ"39  x  6.00  x  IQ23 
T  “f.  x  r2  xTx  10-16 

where  r  is  given  in  Angstroms,  Equation  46  reduces  to 


Bo-  _  67.845 

T  Mj  x  r2xT 

taking  the  log^  one  obtains: 

log  =  1.83152  -log  M|-  log  r 2  -  log  T 
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letting  a  =  2,  substituting  Equation  48  into  Equation  43,  adding  it  to 
Equation  41,  and  dividing  the  result  by  2. 303R,- Equation  49  is  obtained. 

-(9  ~H°-Vot  -t-trons  =  1.5  log  Mi  +  2.5  log  T 
\2.3RT  /  (49) 

—3.5803  +  log  —log  T  +  log  Mj 


Summing  the  electronic  vibrational,  translational,  and  rotational 
contributions  to  the  free  energy  function,  expression  50  is  obtained  for 
the  free  energy  function  of  a  homonuclear  diatomic  molecule. 


(G-H8  Vtrans  +  rot\ 
'2.30RT*  Wvib+elec ' 


=  2.5  log  Mj  +  3.5log  T  +  log 


+  log{Qe|)  -  log  Qv  -  3.5803 


(50) 


In  Equation  50,  Qy  and  are  the  vibrational  and  electronic  partition 
functions  respectively.  Combining  the  above  expression  with  the 
corresponding  expression  for  the  atom  the  change  in  the  free  energy  function 
for  the  reaction  [Ni2  -*■  2Ni]  is  obtained  by  subtracting  the  free  energy 
function  of  Ni2  from  twice  the  free  energy  function  of  Ni  atoms,  which 
gives  the  following  expression  for  the  change  in  the  free  energy  function 
for  the  dissociation  reaction. 


°-s  ,o9  Mi + '-5  iost  -  ,ob 

+  log  Qy  (Nig)  —  2  log  r  (Nig)  +0.3963 


(51) 


Combining  Equation  51  with  Equation  40  and  expressing  the  pressure 
P  in  Torr  one  obtains  Equation  52. 

Jrf  =  -;=3  p^INi>  +  W  +  "’9 


+  iog9li^!ll  +  log  Qv  (Ni2)  +0.3963  +2.8808 
Qe|(Ni2) 


(52) 


where  2.8808  is  the  log  760.  The  conversion  to  Torr  is  required  since 
free  energy  functions  were  defined  with  reference  to  a  standard  state  at 
one  atmosphere. 


51 


AFML-TR-72-217 


The  rotio- 


j%i)  -  P(Ni)Y2 
l+(Ni2)  F(Mi2>Y( 


(53) 


where  y-j  is  the  product  of  the  ionization  cross  section  and  the  electron 
multiplier  efficiency  for  a  given  species.  If  it  is  assumed  as  was  in 
Reference  42  that  (y2/y-|)  =  1,  then  [I+(Ni)/I+(Ni2)]  =  P(Ni)/P(Ni2)  or 
Equation  52  becomes  Equation  54. 


2-4t  =  -  '■">  pT«r«n +  mr + 

+  ,0g  +  ,0g  Qv  (Ni2)+  3.2771 


Utilizing  Equation  54,  the  free  energy  functions  given  by  Hultgren 
(Reference  3)*  the  "best  value"  of  AH2gg  for  the  vaporization  nickel 
obtained  in  this  n=,jort,  the  data  given  by  Kant  for  the  ratio  of  the  ion 
currents,  and  the  constants  in  Table  XIII  a  value  of  Dq(Ni'2)  =  61.7 
±0.6  kcal  per  mole  was  obtained  in  contrast  to  the  value  reported  by 
Kant  (Reference  42)  of  53.3  kcal /mole. 
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APPENDIX  II 

IONIZATION  CROSS  SECTIONS  USED  IN  HASS  SPECTROMETRIC 
VAPOR  PRESSURE  MEASUREMENTS 

Section  1  Examination  of  Ionization  Cross  Section  Data 

In  the  report  of  the  mass  spectrometric  measurements  of  the  vapor 
pressures  and  enthalpies  of  vaporization  of  nickel  oxide  (NiO)g,  nickel 
(Nil  and  oxygen  (09)  from  nickel  oxide  solid  fay  Grimley,  Burns  and 
Inghram (Reference  4)  discussed  on  p.32  etc.,  the  values  of  the  ionization 
cross  section  of  silver  and  oxygen  atoms,  and  the  relative  values  of  the 
ioaization  cross  sections  of  silver,  nickel,  and  oxygen  atoms,  and 
oxygen  molecules  were  required  to  obtain  the  absolute  pressures  of  these 
species  in  the  system.  The  \alues  used  in  Reference  4,  were  the  relative 
cross  sections  calculated  by  Otvos  and  Stevenson  (Reference  51),  which  are 
inconsistent  with  the  experimentally  determined  values  of  cross  sections 
of  silver  and  oxygen  atoms;  therefore,  it  was  necessary  to  re-examine  the 
data  in  Reference  4,  utilizing  the  best  values  of  cross  sections  of 
hydrogen  atoms,  oxygen  atoms,  oxygen  molecules  and  silver  atoms,  and  the 
relative  cross  sections  of  silver  and  nickel  atoms.  The  best  values  of 
these  cross  sections  were  determined  as  indicated  in  the  subsequent 
discussion. 

Since  many  of  the  calculated  and  observed  ionization  cross  sections 
reported  in  the  literature  are  given  relative  to  that  of  the  hydrogen 
atom,  it  was  necessary  to  obtain  a  "best  value"  of  the  cross  section  for 
ionization  of  the  hydrogen  atom.  The  best  value  for  this  cross  section 
was  obtained  from  a  critical  analysis  of  the  ionization  cross  section 
measurements  reported  by  Kieffer  and  Dunn (Reference  52).  Referring  to 
Figure  6  (Figure  3,  Reference  52),  it  is  noted  that  the  two  curves 
representing  ionization  cross  sections  as  a  function  of  the  electron 
energy  are  composites  of  data  of  equal  reliability,  presented  by  four 
different  groups  of  investigators.  The  value  of  the  ionization  cross 
section  which  is  of  primary  interest  in  the  dic'-ussion  is  the  maximum 
value  (crm).  The  "best  value"  of  can  be  obtained  by  averaging  the 
voltage  for  the  maximum  cross  sections,  and  in  turn  averaging  the  cross 
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Figure  6.  Normalized  Cross  Sections  for  Ionization  of  Atomic 
Hydrogen.  O-Fits,  Brackmann  are  the  data  referred 
to  by  the  authors  as  "relative";  they  were  normalized 
to  Born  approximation  calculation  at  500  ev.  X-Fite, 
Brackmann  are  the  data  referred  to  by  the  authors  as 
"absolute";  they  were  normalized  using  the  total  absolute 
cross  sections  measured  by  Tate  and  Sm*th  (Reference  53) 
for  molecular  hydrogen.  References:  Fite,  Rothe, 
Boksenberg,  (References  54,  55  and  56  respectively. 

[Figure  3,  Kiefer  and  Dunn,  Rev.  Mod.  Rhys.  38,  1  (1966)] 
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sections  at  these  voltages.  The  averaging  process  is  valid  since  the 

curves  are  flat  at  the  maxima.  Upon  examinations  one  curve  yielded  a 

-16  2 

value  of  0=  0.66  x  10  cm  at  60  volts  and  the  other  curve  a  value 
«•  -16  2 

of  =  0.70  x  10  cm  at  53  volts,  thus  giving  an  average  value  of 
n*  -16  2 

am  =  0.68  x  10  cm  at  57.5  volts  as  shown  in  Table  XIV  of  this  report. 

The  cross  sections  of  the  oxygen  atoms  and  the  oxygen  molecule  can 
also  be  obtained  from  the  data  available  in  Reference  52.  An  examination 
of  Figure  7  (Figure  6,  Reference  57)  combined  with  the  discussion  on 
page  26,  Reference  52,  leads  one  to  consider  only  the  data  of  References 
55,  56,  and  58.  The  reason  the  data  in  Reference  57  was  not  considered 
was  that  it  was  inconsistent  with  the  data  of  References  55  and  58  which 
were  consistent  with  each  other  in  spite  of  the  fact  that  they  represent 
different  method  of  measuring  cross  sections.  The  data  which  was  considered 

gave  a  value  for  ionization  cross  section  for  0  atom:  o  =  1.56  x  10“16 

2  m 
cnr  at  87.8  volts.  The  ionization  cross  section  for  the  Q2  molecule  as  a 

function  of  energy  is  given  in  Figure  7  (Figure  15,  Reference  52). 

Utilizing  the  discussion  on  page  27,  Reference  52,  a  value  of  o  =  2.80 

-16  2  m 
x  10  cm  at  123  volts  was  obtained  for  02  molecule.  The  above  values 

of  om  for  0  and  02  are  those  shown  in  Table  XIV.  Although  these  values 

are  subject  to  a  number  of  systematic  errors,  there  is  no  evidence  which 

would  lead  us  to  reject  these  results  in  favor  of  any  calculated  relative 

values.  Therefore,  one  is  led  to  use  these  values  as  the  best  available 

ones  when  measuring  pressures  by  mass  spectrometry. 

The  absolute  value  of  the  ionization  cross  section  of  the  silver  atom 

was  measured  by  Crawford (Reference  63)  who  obtained  the  values  shown  in 

Table  XIV.  Although  these  values  disagree  with  those  in  Reference  64 

-16  2 

wherein  the  value  reported  was  a  =  2.68  x  10  cm  at  70  volts,  Crawford's 
value  was  accepted  because  of  the  consistencies  found  between  this  value 
and  other  relative  cross  sections  by  Rovner  and  Norman  (Reference  55). 

For  example,  Rovner  and  Norman  noted  that  if  they  used  their  measured 
value  of  the  relative  ionization  cross  section  of  the  calcium  atom  to  that 
of  the  silver  atom:  [o(Ca)/a(Ag)l  =  0.63,  in  combination  with  absolute 
value  of  the  ionization  cross  section  of  calcium  atom  o(Ca)  =  11.8  x  10"16 
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because  of  mass  analysis  are  cross  sections  tor  tne 
production  of  0+  from  0.  (ira2  =  0.88  x  10"  6  cm  ) 
[Figure  17,  0.  Kieffer  and  G.  H.  Dunn,  Rev.  Mod.  Phys.  38 
1  (1966)] 
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on  obtained  by  McFarland  (Reference  66),  they  obtained  a  value  for 
-16  2 

o(Ag)  =  6.6  x  10  cm  at  50  volts.  Furthermore,  Rovner  and  Norman 

(Reference  65)  noted  that  using  their  measured  value  of  the  relative 

cross  section  of  silver  atom  to  lead  atom,  Ccf(Ag)/cr(Pb)]  =  0.62,  which 

agrees  with  the  value  of  [o(Ag)/o(Pfc)J  =  0.58  reported  In  references 

63  and  67,  and  the  absolute  value  of  o(Pb)  =  8.0  x  10  cm  (Reference  64), 

-16  2 

they  arrived  at  a  value  of  cr(Ag)  =  4.8  x  10  cm  at  50  volts.  Thus  it 

wa*  concluded  in  Reference  65  that  the  beat  value  of  o(Ag)  =  5.6  ±1  x  lO-1'" 

2 

cm  at  50  V.  For  the  purposes  of  this  report  we  intend  to  use  the  best 

-16  2 

available  measured  value  for  Ag,  that  is,  om(Ag)  =  4.63  x  10  cm  at 
70  volts  as  shown  in  Table  XIV,  since  it  is  consistent  with  other  cross 
section  measurements  discussed  above. 

To  obtain  a  relationship  between  the  measured  cross  sections  and 

those  used  in  Reference  4,  we  must  discuss  the  calculation  of  relative 

cross  sections  by  Otvos  and  Stevenson (Reference  51)  and  indicate  the 

limitation  of  this  calculation  in  predicting  ionization  cross  sections. 

Otvos  and  Stevenson  assumed  that  the  maximum  ionization  cross  sections 

are  proportional  to  the  weighted  sum  of  the  cross  sections  of  the  valence 

electrons  of  the  atoms.  The  weighing  factor  was  the  relative  mean 

square  of  the  radii  of  the  electrons.  They  also  assumed  that  the 

molecular  cross  sections  are  the  sum  of  the  atomic  cross  sections.  These 

assumptions  are  limited  by  the  technique  of  calculating  wave  functions 

of  the  atoms  and  also  by  the  fact  that  no  account  is  taken  of  the  energy 

of  the  colliding  electron.  The  effect  of  the  inner  shells  of  the  atom 

on  the  atomic  radii  are  simply  accounted  for  by  correcting  the  atomic 

number  by  a  screening  constant.  The  errors  thus  introduced  become 

accentuated  as  the  atoms  become  more  complex  as  noted  in  a  comparison 

between  the  calculated  relative  cross  sections  of  sodium,  potassium,  and 

hydrogen,  and  the  measured  relative  cross  sections  of  these  atoms.  For 

example,  the  values  of  the  measured  total  cross  section  for  sodium  atom 

at  its  maximum  are:  o_  =  6.8  x  10"^  cm^  (Reference  68);  8.6  x  10"^6 

cm  (References  69  and  70);  and  7.6  x  10  cm  (Reference  71)  which  upon 

-16  2 

averaging  gives  a  value  of  om  =  7.7  x  10  cm  ,  and  a  relative  cross 
section  of  sodium  with  respect  to  hydrogen  (Table  XIV);  [0(Na)/o(H)]  =  11.3. 
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Reference  51  gives  a  calculated  value  of  [o(Na)/o(H)]  =  14.3.  In 

contrast  a  much  larger  discrepancy  between  the  measured  and  calculated 

relative  cross  section  exist  for  potassium.  The  measured  maximum  values 

-16  2 

of  total  cross  section  of  potassium  at  8.5  volts  are  =  7.9  x  10  cm 

-IS  2  m 

(Reference  68),  a  =  9.6  x  10  cm  (References  69  and  70)  and  o_  = 

—16  p  rn  m 

8.6  x  10“  cur  (Reference  71)  yielding  an  average  value  of  crm(K)  =  8.6  x 
-16  2 

10  cm  giving  a  relative  cross  section  with  respect  to  hydrogen  atom 
of  Co(K)/o(H)]  =  12.7,  while  the  calculated  relative  cross  section 
Co(K)/o(H)]  =  38.8  (Reference  51). 

Examining  the  cross  section  for  silver  we  find  similar  discrepancies. 
Utilizing  the  data  in  Table  XIV  it  is  found  that  the  relative  cross 
section  [a(Ag)/cr(H)]  =  6.83  in  contrast  to  the  calculated  value  in 
Reference  51:  [a(Ag)/o(H)]  =  34.8.  Thus  there  may  be  considerable  error 
in  trying  to  use  the  relative  cross  sections  in  Reference  51  to  obtain 
absolute  values.  As  a  comment  on  the  additivity  of  cross  sections  it 
is  noted  that  measured  ratios  of  [a(02)/o(0)j  =  1.80,  which  differs 
from  the  value  of  2.0  predicted  in  Reference  51. 

Since  no  measured  value  of  the  absolute  ionization  cross  section  of 
the  nickel  atom  is  available,  the  relative  value  of  [om(Ag)/om(Ni)] 
reported  by  Cooper  et  al  (Reference  72)  was  utilized  to  estimate  the 
absolute  value  of  om(Ni)  fven  though  it  was  difficult  to  evaluate  the 
experimental  techniques  utilized  by  these  authors.  The  value  of 
[o(Ni)/a(Ag)]  =  0.90  at  SO  V  has  been  reported  in  Reference  72.  Due  to 
the  fact  that  ionization  potentials  for  Ni  atoms  (7.61  V)  (Reference  73) 
are  almost  equal  to  that  of  silver  (7.54  V)  (Reference  73)  it  can  be 
assumed  that  the  ratio  of  ionization  cross  sections  at  the  maximum  are 
approximately  equal  to  the  ratio  at  60  V  or  that  a_(Ni)  =  4.19  x  10"^ 

/>  IN 


The  extrapolation  of  the  ionization  cross  section  as  a  function  of 
energy  is  based  on  the  commonly  used  relationship  between  the  electron 
energy  and  the  ionization  cross  section  where  the  energy  is  less  than 
that  at  which  the  cross  section  is  a  maximum. 
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fv  -  -  (55) 

■'m  Vm-VA 

In  equation  55,  V  is  the  potential  corresponding  to  the  electron  energy, 
and  V  is  the  potential  at  which  the  ionization  cross  section  is  a 
maximum,  while  is  the  appearance  potential,  or  ionization  potential 
for  the  molecular  species  In  question. 

The  above  data  on  cross  sections  can  now  be  applied  to  correct  enthalpy 
and  vapor  pressure  data  reported  by  Grimley,  et  al  (Reference  4)  as  noted 
in  Section  2. 

Section  2  Hass  Spectrometri c  Measurements  of  Vapor  Pressure 

The  measurement  of  vapor  pressure  by  means  of  the  mass  spectrometer 
requires  a  knowledge  of  (1)  the  number  of  particles  of  a  given  molecular 
weight  which  enter  an  ionization  chamber,  (2)  the  fraction  of  particles 
in  the  ionization  chamber  which  are  ionized,  (3)  the  molecular  weight 
distribution  of  the  resulting  ions,  (4)  the  fraction  of  ions  formed  in 
the  ionization  chamber  which  are  accelerated  into  the  magnetic  field  (for 
magnetic  mass  spectrometers)  and  which  arrive  at  the  detector,  (5)  the 
electron  current  provided  by  each  ion  arriving  at  the  detector.  The 
factors  are  often  stated  in  terms  of  a  geometric  or  transmission  factor 
which  defines  the  fraction  of  particles  which  are  transmitted  through 
each  region  of  the  mass  spectrometer;  the  ionization  cross  section  which 
determines  the  fraction  of  particles  that  are  ionized  to  form  a  given  ion; 
and  the  sensitivity  of  the  detector  that  determines  the  detector  current 
produced  per  ion  arriving  at  the  detector.  The  relation  between  the 
ionization  cross  sections,  vapor  pressures,  and  geometric  factors  can  be 
obtained  as  follows.  Referring  to  Figure  8  and  letting 

le  =  electron  current  in  the  ionization  chamber 

A  =  area  swept  out  by  the  electron  beam 

Z  ~  length  of  electron  path  in  ionization  chamber 

o.(V)  =  cross  section  for  ionization  of  species  j  by  the 
J  electrons  of  energy  V 
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Figure  8.  Magnetic  Mass  Spectrometer  Arrangement 


60 


.  .  +. 

(-nr)  =  rate  of  production  of  positive  charge  in  the  ionization 

j  chamber  due  to  species  j 


=  number  density  of  neutral  particles  j  in  the  ionization 
chamber 


The  following  relation  holds  at  constant  electron  energy  between  the 
quantities  defined. 


(56) 


If  the  small  fraction  of  the  molecules  which  are  ionized  can  be 
neglected  the  density  in  the  ionization  chamber  is  obtained  from  the 
following  relation. 

p.  =  flux  out  of  Knudsen  cell  x  geometric  factor  x  (velocity  of 
J  1 

particles) 

p.  =  n.  x  G  x  v-1  (57) 

1 

The  flux  out  of  the  Knudsen  cell  is  given  by  the  following  relation. 
nj  =  Nx  Pj  x  (27rRTM)'1/2  (59) 

where  N  is  Avagadro's  number,  P.  is  the  pressure  of  jth  particle  in  the 

J 

cell,  M.  molecular  weigh.,  T  -  temperature,  R  gas  constant. 

J 

The  molecular  velocity  v  in  the  beam  differs  from  the  average  velocity 
in  the  cell  as  noted  in  Equation  60. 

v  =  (9t7RT/8M),/2  (60) 

The  relation  between  the  mass  flux  from  the  cell,  m,,  and  pressure 

J 

in  the  cell  is  given  by 

1/2 


Pj  =  mj(27rRT/Mj) 


(61) 
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while  the  relation  between  the  n.  and  m.  is  given  by 

J  J 


which  leads  to  the  following  expression  for  p. 

J 

P]  =  GNFj  (2ttRTM)“i/2  (97rRT/8M)',/2 
Rearranging  terms  we  may  write 

G  =  G(27rMr'/2(9;r/8Mr,/2  -(2/3tt)G 
/>  j=  G'NPjtRT)-1 

Tfdnt^  =  G%or)pj  * 

'dt  *  ~  k 

where  k  is  the  Boltzman  constant. 


(62) 

(63) 

(64) 

(65) 

(66) 


To  simplify  the  derivation  of  the  relation  between  electron  multiplier 

current  and  positive  ion  current  consider  the  steady  state  operation  of  a 

continuous  electron  beam  and  particle  flux.  In  this  case,  the  rate  of 

production  of  positive  ions  (dn^/dt)  equals  the  rate  at  which  ions  arrive 

J 

at  the  electron  multiplier  detector  divided  by  a  transmission  factor  (S). 
The  transmission  factor  accounts  for  the  loss  of  ions  between  the  chamber 
and  the  electron  multiplier.  The  electron  multiplier  current  Ip  is 
proportional  to  the  ion  current  arriving  at  the  multiplier  and  is  given 
by  Equation  67. 

'p  =  T,  (f)?  =  FTi’i  1671 

F  is  the  electron  multiplier  amplification  factor,  and  y^  is  the 
sensitivity  of  the  electron  multiplier  for  the  ion  j.  The  sensitivity  is 
primarily  a  function  of  molecular  weight  ,  for  ions  of  atoms  and  those 

o 

of  simple  molecules  of  the  same  energy.  If  the  simple  ion  fragments  at 
the  electron  multiplier  surface,  the  sensitivity  is  then  a  function  of  the 
momentum  of  the  fragments. 
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Equation  66  can  now  be  rewritten  as: 

i+t  -  ML  -  G  le^°‘i  *pi  S* 

1  j  FYj  "  k 

where  k  is  the  Boltzman  constant. 


(68) 


Rewriting  Equation  68  by  taking  logarithms  to  the  base  e,  (In), 
Equation  69  is  obtained: 


Inlj+T=  In  le  +  Ino-j  +  InPj  +  ln(G*Sf/k) 

(69) 

In  lj"T  —  In  <r  j  =  In  Const  +  In  Pj 

(?Oo) 

In  Ipj  T  -  Incrj  =  In  Const*  +  In  Pj 

(70b) 

It  is  seen  from  the  above  relation  that  errors  in  the  absolute  cross 
section  directly  affect  the  value  of  the  (Const)  and  subsequently  all 
values  of  the  measured  pressures  which  are  deduced  by  using  the 
calibration  constants.  Since  Equation  55  is  used  to  obtain  c,  from 
[ajm]  the  errors  in  the  use  of  the  maximum  cross  section  given  by  Otvos 
and  Steverson  (Reference  51)  reflect  directly  on  the  value  of  a-  used 

J 

in  Equations  70a  and  70b,  and  subsequently  on  the  deduced  absolute 
pressures.  For  example,  the  value  of  the  cross  section  of  silver  which 
is  deduced  on  the  basis  on  the  Otvos  and  Stevenson  calculation  is 

a  9  I  r  o 

34.8Tra4  where  Traf  =  0.879  x  10"  0  cm4  is  the  cross  section  of  the  H  atom 

U  U  -16  ? 

or  am( Ag,  70  ev)  =  30.6  x  10  cm.  The  measured  value  is  of Ag)  = 

4.63  x  10  cm4.  See  Table  XIV.  The  calibration  data  given  in 

Reference  4  indicates  that  constant  in  Equation  70a  in  error  by  a  factor 

of  -ln(30.64/4.63).  Furthermore  the  relative  cross  section  of  Ni  is  given 

in  References  4  and  51  as  24.4  which  gives  an  absolute  value  of  a_(Ni)  == 
-16  ?  ' 

21.5  x  10  cm  in  contrast  to  the  value  obtained  from  the  combination 

1  c  p 

of  several  measurements,  i.e.  [om(Ni)]  =  4.19  x  10  cm  .  The  error 
in  the  cross  section  of  nickel  atom  partially  compensates  for  the  error 
in  the  (constant)  in  Equation  70b.  For  example,  the  relation  between 
the  constants  reflecting  the  use  of  two  different  cross  sections  for 
silver  is  given  in  Equation  71. 
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In  [(const  Jj/tconstg)]  =  In  (o^/a-j^g  (71) 

where  const-j  is  the  constant  obtained  in  Reference  4. 

Similarly  utilizing  the  different  cross  sections  for  nickel,  one 
obtains  from  Equation  71  an  expression  for  correcting  the  data  of 
Reference  4  by  the  use  of  the  measured  cross  section.  Rewriting 
Equation 70a  to  apply  to  the  nickel  ion  current,  Equation  72  is  obtained: 

In  TI+(Ni)-ln  cr(Ni)  =  In(const)  +  InP(Ni)  (72) 


If  we  designate  the  vapor  pressure  measured  in  Reference  4  as 
[P(Ni)].j  and  the  corrected  vapor  pressure  as  [P(Ni)]2  Equation  73  is 
obtained. 

InJcrgfNiJ/cj-,  Wo]  =  In  (const  | /const 2)  +  In[P|(Ni)/P2(Ni)j  (73) 

Substituting  ln[a2(Ag)/a-j(Ag)]  for  ln[const^)/(const2)]  from  Equation  71, 
Equation  74  is  obtained  which  leads  to  a  relation  between  the  reported 
and  the  corrected  vapor  pressure.  Equation  75. 

ln[a2(N«)/o-,(Ni)]  =  ln^(Ag)/o-|(Ag)]  +  ln[p,(Ni)/P2(Ni)]  (74) 

ln[o-2(Ni)/o-|(Ni)][CT1(Ag)/<72(Ag)]  =  ln[p,(Ni)/P2(Ni)]  (75) 

Utilizing  the  values  for  given  in  Table  XIV  the  corrected  value 
P2(Ni)  =  0.78  P^(Ni)  is  obtained. 

Grimley,  et  al  (Reference  4),  also  reported  the  vapor  pressure  of 
oxygen  in  equilibrium  with  NiO.  This  data  was  examined  in  light  of  the 
absolute  cross  section  data  that  is  available.  They  determined  the 
pressure  of  oxygen  from  a  measurement  of  the  ratio  of  0+  ion  to  the  02 
which  was  assumed  to  be  related  by  the  equilibrium  constant.  Equation  77 
for  the  following  reaction: 


I/2C>2  0 

(76) 

K  = _ EIQ1 _ 

P  [P(02)]'/2 

(77) 
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A  value  for  the  pressure  inside  the  Knudsen  cell  can  be  obtained 
from  the  following  consideration;  assuming  that  the  reaction  in  the  cell 
is  given  by  Equation  78  and  that 

NiO^  -  =7  Ni(s)  +  l/20g  (78) 

the  pressure  of  oxygen  in  the  cell  must  satisfy  the  following  steady 
state  relationship:  [Rate  of  Vaporization  of  Og]  =  [Rate  of  Condensation 
of  Og]  -  [Rate  of  Effusion  of  0^]  -  [1/2  Rate  of  Effusion  of  0].  (79) 

The  last  term  accounts  for  the  loss  of  0g  due  to  the  dissociation 
Reaction  76  and  can  be  related  the  partial  pressure  of  0  atom.  If  the 
rate  of  vaporization  is  the  maximum  possible  for  the  temperature  of  the 
system  then  it  is  equal  to  the  maximum  rate  of  condensation,  that  is  the 
rate  of  collision  of  02  molecules  with  the  sample  surface  at  the 
equilibrium  pressure. 

The  individual  rates  noted  above  can  be  expressed  as  follows  if  we 
assume  that  the  vaporization  coefficient  and  Clausing  factors  are  unity 

Rote  of  Vaporization!  =  A  [p(02)]e[25rkTM(02) /n]"I/2  (80) 

Rote  of  Condensation:  [— ^•]coa  =  A[p(02)]c  [2rrkTM(02)/N]“l/2  (81) 

Rote  of  Effusion  of  02:  =  a[p(02)]e [27rkTM{02)/fsl]— 1/2  (82) 

Rate  of  Effusion  of  O:  [^5r^]eff  =  a[p(0)]c[27rkTM(0)/N]_l/?  (83) 
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In  the  above  expression  A  is  the  cross  sectional  area  of  the  sample 
and  a  is  the  cross-sectional  area  of  the  effusion  hole.  For  any  given 
species  the  relation  between  the  equilibrium  pressure  and  the  pressure 
in  the  cell  is  obtained  from  the  above  relations  which  reduce  to 

A[p{02)]e  =  a[p(02)]c  +  s[p(°2^c  +  I/2o[p(0)]c[m(02)/M(°)]I/2  (84) 

The  ratio  [M(02)/M(0)]^Z  =  /T while  the  ratio  [P(0)]/[P(02)]^Z  =  K. 
Noting  that  (a/A)  &  .01  and  the  P(0)£J  .01  P(02)  we  may  neglect  the 
difference  between  [P(02)]c  and  [P(02)]g  and  further  assume  that  mass 
flow  rates  from  the  cell  are  equal  to  those  at  the  equilibrium  pressures. 
These  relations  allowed  Gri ml ey.et  al  (Reference  4)  to  use  the  observed 
oxygen  ion  currents  to  determine  the  constants  in  equation  70,  by  using 
the  ratio  of  cross  sections  as  noted. 


ln[TI+(0)/TI  +  (02)]  -ln[<r(0)/<r(02)]  =  In [p(0)/P(02)] 

[$§M^HPI0)/PI02>] 


(85) 

(86) 


Substituting  the  expression  for  the  equilibrium  constant  for  the 
ratio  of  the  pressures, 

[l+(0)/l+(02)J[cr(02)/a(0)]  =k[p(02)]  1/2  (87) 

The  ratio  [a(Q)/c(Q2)3  =  0.50  if  the  cross  sections  in  Reference  51 
are  used  and  the  [a(0)/o(02)]  =  0.56  if  the  data  in  Table  XIV  is  used. 
The  value  of  the  ratio  [P-| (02)/P2(02)]  which  is  determined  by  the  use 
of  the  Otvos  and  Stevenson  (Reference  51)  cross  sections,  in  conjunction 
with  the  measured  cross  sections  is  given  by 

(02)/o'|(0)j/^or2(0)/cr2(02)j=  jV2(02)/P|(02)J  (88) 

or  that  the  pressure  of  oxygen  as  determined  by  using  the  measured 
instead  of  the  calculated  cross  section  is  1.24  [P(02)]  as  determined  by 
the  use  of  the  calculated  (0  +  S)  (Reference  51)  cross  section]. 
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Since  the  (const^)  in  equation  70a  assumes  that  the  electron  beam 
current  is  a.  constant,  no  relation .between  the  apparatus  (constant)  as 
determined  by  the  oxygen  calibration  technique  vs.  the  silver  calibration 
can  be  readily  drawn.  However,  it  is  noted  that  under  the  same 
experimental  conditions  the  use  of  the  calculated  cross  sections  instead 
of  the  measured  cross  sections  introduces  a  ratio  of  3.2  between  the 
factors  in  the  (constant)  which  are  independent  of  the  operating 
conditions  of  the  instrument. 

Since  the  pressure  of  NiO^  is  of  interest  in  determining  the  heat 
of  vaporization  of  Ni0^s ^ ,  the  pressures  reported  by  Grimley  (Reference  4) 
must  be  corrected  for  the  change  in  cross  sections  reported  in  Table  XIV. 
Grimley,  et  al,  reported  using  a  relative  cross  section  27.7  for  NiO  which 
is  equal  to  the  sum  of  the  (0+S)  cross  sections:  o(Ni)  +  o(0).  Due  to 
the  lack  of  more  exact  information  in  this  report  it  is  assumed  that 
a(NiO)  =  a(Ni)  +  a(0),  but  that  the  value  of  cr(Ni)  =  4.19  x  10~16  cm2 
and  o(0)  =  1.56  x  10~16  cm2  or  the  value  cr(NiO)  =  5.75  x  10"16  cm2.  The 
correction  to  the  pressure  of  (NiO)g  due  to  the  correction  in  cross  section 
is  given  by  Equation  75  as  applied  to  NiO  where  o1  is  the  (0  +  S) 

(Reference  51)  cross  section  and  o2  are  the  measured  cross  sections  in 
Table  XIV. 

r^WiOh  rcr,  IAg)n  _  rP^NiOh 

Lct|  (NiO)  J  Lo-2(Ag)J  ‘  LP2(NiO)  J 

Substituting  the  values  of  the  cross  sections  from  Table  XIV  into 
Equation  89  gives  P2  =  0.641P^.  This  value  of  P2  was  used  to  recalculate 
the  heats  of  vaporization  of  NiO^  as  reported  in  Part  II  of  this  report. 
Corrections  were  also  made  to  reported  pressures  of  02  and  the  data  of 
Reference  4  was  recalculated  (Part  II  of  this  report).  The  recalculation 
indicated  that  the  mass  spectrometer  is  reliable  under  suitable  operating 
conditions  in  generating  absolute  pressures. 
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TABLE  XIV 

IONIZATION  CROSS  SECTIONS 


Species 

Voltage 

0  x  10+16  cm2 

°m(0tS)f 

Exp  relative 

H 

Max  57.5(a) 

0.68^ 

1 

1 

0 

Max  87.8^ 

1.56{b) 

3.29 

2.29 

°2 

Max  123.0(^ 

2.80^ 

6.58 

4.12 

Ag 

Max  72^ 

4.63<d> 

34.8 

6.83 

Ag 

50 

4.48(a) (b) (c) (d) (e) (f) 

Ag 

60 

4.50<d> 

Ni 

60 

4.10(e) 

Ni 

Max  70 

4.19 

24.4 

(a)  Figure  7  (Reference  52,  p.  11,  Figure  3). 

(b)  Figure  8  (Reference  52,  p.  12,  Figure  6,  data  of  Rothe  and  Fite). 

(c)  Reference  52,  p.  18,  Figure  17.  The  total  ionization  cross 

section  was  used  for  these  calculations  because  the  electron 
energy  used  in  mass  spectroi.'etric  studies  may  vary  from  10  to 
70  ev  but  it  is  often  limited  to  20  ev  when  dissocation  is  to 
be  avoided.  In  cases  where  the  electron  energy  is  greater  than 
20  ev  corrections  should  be  made  for  the  dissociation  of  0?. 
Reference  52,  p.  43.  c 

(d)  Reference  63,  Figure  2. 

(e)  Combination  of  data  References  63,  65,  67,  and  72. 

(f)  Reference  51. 
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SUMMARY 

The  Langmuir  technique  for  determining  vapor  pressures  was  utilTZecI 
to  determine  the  vapor  pressure  of  nickel  between  1233  and  1683°K.  The 
vapor  pressure  data  thus  obtained  was  combined  with  other  available  data 
in  the  literature  to  obtain  a  best  value  for  the  heat  of  vaporization  of 
nickel  solid  and  liquid  between  the  temperatures  1233°K  and  1895°K. 
Equations  for  the  vapor  pressure  of  nickel  in  these  temperature  ranges 
were  also  determined. 

The  data  in  the  literature  used  to  determine  heats  of  dissociation 
of  the  NiO  and  Hi 2  gaseous  molecules  was  corrected,  and  revised  heats  of 
dissociation  of  Nig  and  NiO  gaseous  molecules  were  obtained,  as  well  as 
estimates  for  heat  of  vaporization  of  NiO  solid. 

The  literature  on  ionization  cross  sections  useful  in  mass  spectrometry 
was  examined,  and  best  values  of  the  ionization  cross  sections  of  silver, 
lead,  and  nickel  were  determined. 
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